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Abstract 

Given a finitely generated semigroup S of the (normed) set of lin- 
ear maps of a vector space V into itself, we find sufficient conditions 
for the exponential growth of the number N(k) of elements of the 
semigroup contained in the sphere of radius k as k — > oo. We relate 
the growth rate lim^oo log N(k) j log k to the exponent of a zeta func- 
tion naturally defined on the semigroup and, in case S is a semigroup 
of volume-preserving automorpisms, to the Hausdorff and box dimen- 
sions of the limit set of the induced semigroup of automorphisms on 
the corresponding projective space. 

1 Introduction 

The asymptotic behaviour of the norms of products of some fixed finite set of 
square matrices has been extensively studied in the context of the theory of 
random matrices. In particular, in a celebrated paper [FK60], Furstenberg 
and Kesten proved that, given some finite number of square matrices A t , 
under suitable conditions the norm of almost all products of k of the Ai 
grows as r y k , where 7 is the Lyapunov exponent associated to the A^. In 
this paper we address the subject from a different point of view, namely we 
consider all possible products of the Ai and provide sufficient conditions for 
the existence and boundedness of Hindoo log N(k)/ log k, where N(k) is the 
number of these products that lie inside the (closed) sphere of radius k. As 



Figure 1: (left) The Cubic gasket C3 C M.P 2 in the triangle T of vertices with 
homogeneous coordinates [x : y : z] = [0 : : 1], [1 : : 1], [0 : 1 : 1]. The 
picture shows (in green) the set T7C3 i n the affine chart z = 1. (right, top) Log- 
log plot of the norms of matrices Cj £ C3, |/| < 11, ordered in lexicographic 
order. The fastest growing norms are ||Cj ■ Cj+i • • • Cj+fc|| ~ a*, where sums of 
indices are intended "modulo 3" and ~ 1.84 is the Tribonacci constant. The 
slowest growing ones are = k. (right, bottom) Log- log plot of the function 

Nc 3 (k) representing the number of matrices of C3 whose norm is not larger than 
k. Numerical data (the values of Nc 3 (k) shown in the graphic are exact, see 
Table 3) indicate that N Ca (k) ~ Ak s for A ~ 0.967 and s ~ 2.444. According to 
Conjecture 1, this entails that dime C3 > 1.63. 

a byproduct, we relate the rate of this growth to a zeta function naturally 
denned on semigroups of square matrices and, in particular cases, to the 
Hausdorff dimension of the limit set of the orbit of a point under the natural 
action induced by the semigroup on its corresponding projective space. 

Motivational Example 1: The Cubic gasket C 3. The real self-pro jective 
fractal C 3 C RP 2 (see Fig. 1) is the main reason for our interest in the sub- 
ject of the present paper. It was first introduced, in the author's knowledge, 
by G. Levitt [Lev93] and independently rediscovered more recently by the au- 
thor and LA. Dynnikov in connection with the S.P. Novikov theory of plane 
sections of periodic surfaces [DD09] . We call it cubic because it is related to 
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the topology of plane sections of the cubic polyhedron {4,6|4} (see [DD09] 
for details) and gasket because it has the same topology of the Sierpinski and 
Apollonian gaskets. Like the Sierpinski gasket, it can be thought as the set 
obtained by removing from the (projective) triangle T(£) with vertices [ei], 
[ e 2], [^3], where £ = {ei,e 2 ,e 3 } is any frame of M 3 , the (projective) trian- 
gle with vertices [ex + e 2 ], [e 2 + e 3 ], [e 3 + ei] and repeating this procedure 
recursively on the three triangles left. We denote by Tf~ : c 3 C T(£) the set 
obtained after repeating this procedure k times. Clearly C3 = H^T^c^, 
i.e. we can get as close as we please to C3 by considering sets Tk t c 3 with 
large values of k. C 3 can also be characterized as the (unique) subset of the 
triangle with vertices [1 : : 0], [0 : 1 : 0] and [0:0:1] which is invari- 
ant under the action of the (free) subsemigroup of PSL 3 (N) generated by 
the projective automorphsims fy, i = 1,2,3, induced by the following three 
SL 3 (N) matrices: 

(\ o\ /11 o\ A i\ 
d = 1 1 , c 2 = 1 , c 3 = 1 1 . 
\i 1/ \o 1 1) \o 1/ 

By abuse of notation, we denote by C 3 also the semigroup generated by the 
Cj. As a consequence of a conjecture of S.P. Novikov [NovOO], the set C3 
is supposed to have Hausdorff dimension strictly between 1 and 2. What 
makes checking this conjecture non-trivial is that each ipi has exactly one of 
the three vertices as fixed point and in that point it has Jacobian equal to 
1 3 , namely the iterated function system (IFS) {^i, ^2, ^3} is parabolic rather 
than hyperbolic 1 . 

No analytical bound for this fractal is known to date. In Section 4, 
based on numerical evidence and Theorems 4 and 5, valid for 2 x 2 matrices, 
we conjecture that the box dimension of C 3 is related by the growth rate 
of the norms of the elements of the semigroup generated by the (see 
Fig. 1), namely by s = lim^oo log N(k)/ log /c, where N(k) is the number of 
matrices of C3 inside the closed ball of M 3 (R) of radius k. According to this 
conjecture, dims C3 > 2s/3 ~ 1.63 (see Section 4.3). 

Motivational Example 2: The Apollonian gasket. The complex self- 
projective fractal A3 C CP 1 is possibly the fractal with the oldest ancestry, 

1 Recall that a IFS {/1, . . . , f m } on a metric space (M, d) is said hyperbolic when all /, 
are contractions with respect to d and parabolic when all /j are non-expanding maps. 
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since its construction relies on a celebrated result of the Hellenistic mathe- 
matician Apollonius of Perga (ca 262 BC - ca 190 BC), known in his times 
as The Great Geometer. Apollonius' result, contained in the now-lost book 
Tangencies but fortunately reported by Pappus of Alexandria in his Col- 
lection [Pap40] published about five centuries later, concerns the existence 
of circles tangent to a given triple of objects that can be any combination 
of points, straight lines and circles. In particular, given three circles which 
are mutually externally tangent to each other (sometimes called the four 
coins problem [01d96]), there exist exactly two new circles tangent to each 
of them, one externally and one internally (see Fig. 2). The three given 
circles plus any one of the new ones 2 form a Descartes configuration, since 
it was Descartes that stated the following remarkable relation between the 
curvatures ci,...,C4 of the four circles (see [Cox37] for details) memorial- 
ized three centuries later by the Chemistry Nobelist Frederick Soddy in 
his poem "The Kiss Precise" [Sod36] after rediscovering it independently: 

2£ti c ? = (£li c *) 2 - 

Since Mobius transformations preserve circles and are transitive on triples 
of distinct points, they also act transitively on the set of all possible Descartes 
configurations; this fact suggests that their most natural environment is the 
Riemann sphere CP 1 rather than the plane. Any Descartes configuration D 
divides CP 1 in 4 curvilinear triangles in such a way that every circle of 
D is one of the two Soddy circles of the remaining three circles of D. By 
drawing the new Soddy circle of each of the 4 triples we are left with 4 new 
Descartes configurations. By repeating this process recursively we generate 
an infinite osculating circle packing of CP 1 which, not surprisingly, is called 
Apollonian packing. 

Here we rather focus our attention on any one of the curvilinear triangles 
T and call Apollonian gasket A3 the set of points of T left after removing from 
T the interior of all Soddy circles inside it. Like in case of the cubic gasket, 
A 3 can be characterized as the invariant set of a complex self-pro jective 
parabolic IFS. The fact that, thanks to the complex structure of CP 1 , A 3 is 
self-conformal was exploited by Mauldin and Urbanski to prove some of its 
fundamental properties [MU98] . Unfortunately these techniques do not seem 
to extend to the previous (real) case, when the IFS maps are parabolic but 
not conformal. 

In 1967 K.E. Hirst [Hir67] introduced the Hirst semigroup H, namely the 
2 In Soddy 's honor the two new circles are called Soddy 's circles. 
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subsemigroup of SX 4 (N) generated by the matrices 
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as an effective tool to represent the radii of the Soddy circles in the gasket. 

In a series of fundamental contributions to the study of the Hausdorff 
dimension of the gasket [Boy 70, Boy71, Boy72, Boy 73a, Boy 73b, Boy82], 
D.W. Boyd ultimately characterized this dimension in terms of the Hirst 
semigroup by proving (implicitly, in terms of the circles' curvatures) that: 1) 
the number Nn(k) of the semigroup matrices with norm 3 not larger than k is 
logarithmically asymptotic to k s for some s > 0, namely lim^oo log 1 ^ ( - fc - > = s; 
2) diniH- A 3 = s. 

Later in this paper we show that A3 can be seen as the invariant set 
of the parabolic Kleinian IFS corresponding to the subsemigroup of ST^C) 
generated by the matrices 

which by abuse of notation we will denote too by A3, and conjecture that 
linifc^oo log Na 3 (k)/ log k = 2 dim// A3 based on the fact that this relation 
holds for similar semigroups that induce hyperbolic IFSs and on numerical 
evidence. 

Motivational Example 3: The Sierpinski gasket. The self-affine fractal 
set S3 C M 2 is less ancient than the Apollonian one, having been introduced 
in the Mathematics literature by W. Sierpinski only in 1915 [Siel5], but it 
does have a long history too since its pattern has been known and used in art 
for about a millennium [PA02] (see Fig. 3). Its dimension is easily calculated: 
dim H S 3 = log 2 3 (e.g. see [Fal90]). Since both PSL 3 (R) and PSL 2 {£) 
contain a subgroup isomorphic to the group of affine transformations of the 
plane, the Sierpinski gasket can also be seen as a real (respectively complex) 
self-projective fractal of IRP 2 (respectively CP 1 ). 

3 Since all norms are equivalent in finite dimension, this is true for any norm. 
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Figure 2: (top, left) Inscribed and circumscribed circles in the four coins problem, (bot- 
tom, left) Apollonian gasket A3 in the curvilinear triangle with mutually tangent arcs as 
sides and the points of homogeneous coordinates \z : w] equal to [1 : 1], [—1 : 1], [i, 1] as 
vertices, represented in the affine chart w = 1. In the image it is shown (in green) the 
set TV Ag. (right) Log-log plots of the norms of the matrices of the semigroup H in lexi- 
cographic order (top) and of the relative function Njj(k) (bottom) counting the number 
of matrices of H whose norm is not larger than k (see Table 3 for the values of Nnik) 
shown in the graph). 



6 




Figure 3: (left) Image of the Sierpinski gasket in the triangle T with vertices (0, 0), 
(1,0), (0, 1). In the picture it is shown, in green, the set 7Vs 3 - (right) Detail of 
a cosmatesque [PA02] mosaic dated about llth-12th century (photo taken by the 
author at the Phillips Museum in Washington DC). 



A semigroup generating the Sierpinski fractal in MP 2 is, for example, the 
one generated by the matrices 




Sf = ^= 1 1 ,5» = -^= 2 0\,S : 





'3 



One generating it in CP 1 is, for example, the one induced by the matrices 

c _ Wl A c 1/1 l\ c 1/1 -1 
5l " J2 V0 V ' 2 " V2 \0 2) ' ^ 3 " V0 2 

The regularity of the matrices S* and Sf makes possible to perform simple 
direct calculations that illustrate the main points of this paper. 

Consider first the real version. Let H^Hoo be the norm given by the 
maximum absolute row sum of S. Since all lines of the Sf~ sum to 2 then 
||5« • • • SfJ\ = (2" 1 / 3 ) p • 2P = 2 2 p/ 3 for every p > 1. Hence in the sphere of 
radius k lie 

L| ^ff J 3 L|io g2fc +ij _ x 
N(k)= 3P = 2 

p=0 
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products of the Sf, where |_31og 2 /c/2j is the integer part of 31og 2 /c/2, and 
therefore 

s R = hm — — = - log 2 3. 

fc^oo log 2 k 2 

Note that sr is also the exponent that separates the values for which the 
series 2~2se(s R ) \\^\\~ s diverges from those for which it diverges, where the 
sum is extended to all elements of the semigroup freely generated by the 
Sf. Finally, note that the following relation holds between the Hausdorff 
dimension of the Sierpinski gasket and the rate growth: 3dim# S 3 = 2s R . 

Consider now the complex version. Endow M 2 (C) with the norm \\S\\ 
given by the largest modulus of the entries of S. Since the last row of each 
Sf is (0, 2), then ||S£ • • • Sf p \\ = (2~ 1 / 2 ) p • 2 P = 2 p / 2 for every p > 1. Hence in 
this case 

[2log,k\ oL21og 2 fe+lj _ i 

N(k)= £ 3P = " 2 

and therefore 

sc = hm — — = 2 log 2 3. 

fc^oo log 2 k 

Similarly to what happens in the real case, sq is also the exponent that sep- 
arates the values for which the series ?~2se(s c ) \\^\\~ s diverges from those for 
which it diverges, where the sum is extended to all elements of the semigroup 
freely generated by the Sf. Note that in this case the relation between the 
Hausdorff dimension of the Sierpinski gasket and the norms' growth rate is 
the following: 2 dim// S3 = Sc- 

For thorough surveys on the Sierpinski gasket and especially on the 
more challenging Apollonian gasket we refer the reader to the book by 
A. A. Kirillov [Kir07], the series of papers by Lagarias, Mallows, Wilks and 
Yan [GLM+03, GLM+05, GLM+06] and the recent article by Sarnak [Sarll]. 

The present paper is structured in the following way. 

In Section 3 we generalize Boyd's arguments on the asymptotics of the 
sequence of radii of Soddy's circles in a Apollonian gasket and use them to 
obtain similar results on the asymptotics of norms of matrices in subsemi- 
groups S C M n (K), K — R, C, by introducing a sufficient condition for the 
existence of lim^oo log Ns(k)/ log /c, where N$(k) is the number of matrices 
in S whose norm is not larger than k, and relating this limit to the critical 
exponent of a natural zeta-function defined on S. 
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In Section 4 we consider the action induced by these semigroups on the 
corresponding real or complex projective spaces and study, in particular but 
significant cases, the relation, observed above in case of the Sierpinski gasket, 
between the critical exponent of the semigroup and the Hausdorff (for n — 2) 
or box dimension (for n > 3) of the limit set of a point under its action. 

We use Section 2 below to introduce the main concepts, notations and 
definition used throughout the paper and to state the main results of the 
paper. 



2 Notations, definitions and main results. 

Matrices and Norms. We endow the vector space M n (K) of all n x n 
matrices with coefficients in K with the max norm, namely, given a matrix 

M = (Mj), 

We denote by B™ C M n (K) the closed ball of radius r > in this norm. 
Note that this norm is not sub-multiplicative but rather 



SU P npim^n = n - (*) 



JPQW 

P,QeM n (K) 

Since in finite dimension all norms are equivalent, the main results of the 
paper will not depend on this particular choice. 

The multi-indices semigroup. We denote by X m the infinite m-ary tree 
of multi-indices of integers ranging from 1 to m defined as follow. The root 
of the tree is the number 0. The to children (1-indices) of are the integers 
from 1 to to. Their children (2-indices) are the ordered pairs li, . . . ,mi and 
so on recursively for the k-indices, k > 2. We denote by X™ the set of all 
/c-indices of X m . Since we will use them often, we denote by T)™, £ > 0, 
the set of all diagonal multi-indices I — %\. . .% k G X m , k < £, i.e. such that 
%x — ■ • • — i k , and set V m = U£> X > ™. Similarly, we denote by J e m , £ > 2, the 
set of all next-to- diagonal multi-indices / = . . . i k e X m , k < £, i.e. those 
such that %\ 7^ %i — ■ • • — ifc, and set J m = U^> 2t 7/ n - 

We endow X m with the canonical structure of semigroup given by i± . . . % k ■ 
i\ . . . i' k , — %\ . . . i^i! x . . . i' k , with as identity element. We also endow X m with 
a partial order by saying that I > J if I can be factorized as / = LJ for some 
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multi-index L ^ 0. Finally, we denote by I' — i\ . . . i& the /c-index obtained 
from the (k + l)-index I = i ii . . . by dropping the first index on the left. 



Gaskets of matrices. Given m matrices A 1 , . . . ,A m e M n (K), K = R or 
C, we denote by A = (Ai, . . . , A m ) the semigroup they generate, given by 
the intersection of all subsemigroups of M n (K) containing all the Ai and the 
unit matrix. 

In this paper we are mainly interested in the asymptotic growth of norms 
of matrices in free semigroups but, since our results hold for the more gen- 
eral case when there are relations between the generators, we often formu- 
late theorems using the more general concept of semigroup homomorphisms 
X m — > M n (K). We often denote such objects with the letter A and use the 
notation 



where Ai = A{i). We say that the matrices A ± , . . . , A m generate A. 

Notice that, when there is no relation between the A i: then there is a 
bijection between A(X m ) and {!„} U A, so when the Ai are free generators 
it is essentially equivalent referring to either the semigroup homomorphism 
A or the semigroup A. 

Given any M e GL n (K) we denote by Am the "right coset" map defined 

by Am (I) ==A/M. If A is free then there is a bijection between Am^ 171 ) and 
{M} U AM, where AM is a right coset of A. Clearly A ln = A. 

Definition 1. We denote by N^ M (r) the cardinality of the set B^nAM(1 m ) 
and say that Am is a m-gasket ( or simply a gasket ) if N^ M (r) < oo for 
every r > 0. We say that the gasket Am is hyperbolic if the sequence = 
min/exm ||A/M || diverges exponentially, namely if there exists a > 1 such 
that Ok x a k , where x means that the ratio of the terms on either side is 
bounded away from and oo for all k. When is slower than exponential 
we say that Am is parabolic 4 . 

Example 1. Every semigroup A : X m — > GL n (C) whose generators have all 
their spectrum outside the unit circle is a hyperbolic gasket. Consider for 
example the simple case of A : X 2 — >■ GL 2 (C) with 



4 Note that au cannot be faster than exponential so this covers all possible cases. 



A I ^A(i 1 ...i k )=A 
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where |A| > 1. Then mm^^k ||4r|| = ||A^|| = ||A 



Example 2. Every semigroup A : X m — > GL n (C) whose generators have all 
non-zero coefficients strictly larger than 1 is a hyperbolic gasket, since the 
norm of k of such matrices will be not smaller than the k-th power of their 
smallest non-zero entry. 

Example 3. The (free) semigroup C 2 C SL 2 (N) generated by the two parabolic 
(i.e. with trace equal to ±2) matrices 



is a parabolic gasket. It is a gasket because, if M e SL 2 (N) is distinct from 
the identity, then \\CiM\\ > ||M|| + 1 since M has at least a column with two 
entries different from 0. It is parabolic because 



Note that C 2 contains also hyperbolic elements (namely matrices Cj with 
\trd\ > 2), e.g. 



so that in the sets C 2 ^ = {Ci, |/| = k} there are some elements whose norm 
grows polynomially and some others whose norm grows exponentially with k, 
similarly to what happens for the cubic and Apollonian gaskets (see Figs. 1 
and 2). 

This elementary but still non-trivial example was suggested to the author 
by LA. Dynnikov and was the starting point for the author's study of the 
asymptotics of norm's growth in semigroups of linear maps in full generality. 

Example 4. Suppose that Ay, . . . , A m e M n (Z), where Z = Z or Z[i], gen- 
erate freely A. Then A is a gasket, since in M n (Z) there are only finitely 
many matrices whose norm is smaller than any fixed r > and by the freedom 
hypothesis the products of any number of Ai are all distinct. 

Note that clearly if A is a gasket then so is also every semigroup conju- 
gated to it, as well as every semigroup obtained from it by multiplying all 
elements by a constant A such that |A| > 1. 




mm{||C7||}<||^ 



= k. 




11 



Definition 2. Let A be a n-gasket and M E GL n (K). We call zeta function 
of Am the series 

Ca m (s) = \\A r M\\ s ' 
We call exponent of Am the number s_a m defined as follows: 

SA M = SUp{s|C4 M ( s ) = °°}- 

s>0 

Note that, if sa m < oo, we also have that sa m = inf s > {s | Ca m ( s ) < °°}- 
Example 5. Let A be the (parabolic) gasket generated by 

A x = j), A 2 = A 1 eSL 2 (N). 

In this case \\A Ik \\ = k, so that Ca( s ) — J2k&N^ k ^~ s diverges for all s, i.e. 
Sj[ = oo. 

On the contrary, let B be the (hyperbolic) gasket generated by 

B v = l ^j,B 2 = B 1 e GL 2 (N). 

Then \\Bi k || = F 2 k+ 2 where F = (0, 1, 1, 2, 3, 5, ... ) is the Fibonacci sequence. 
Hence asymptotically \\Bj k \\ ~ g 2k , where g = 1+ 2 V ^ is the golden ratio, and 
so Cb(s) diverges or converges with ^2 ken 2 k g ~ 2sk , i-e. S& = 2 lo * - . 

Next proposition shows that norms of matrices in Am have the same 
asymptotic properties of those in A for every M e GL n (K): 

Proposition 1. Let A be a m-gasket and M G GL n (K). Then Am is a 
m-gasket and s^ M = s^. 

Proof. It is a direct consequence of n\\P\\ ||M|| > \\PM\\ > n ff-i\\ - D 

Showing that hyperbolic gaskets have a finite exponent does not require 
any effort: 

Proposition 2. Let A : T m — > M n (K) be a hyperbolic gasket. Then is 
finite. 
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Proof. Since A is hyperbolic then \\Aj\\ > Act^ for some A > and a > 1. 
Hence 

oo oo 

Ca(s) < J^m^'a* 8 = ^A s m fc(1 - slog ™ a) , 

fc=0 fc=0 

so that (a{ s ) < oo for s > log Q m, namely < log a m. □ 

As Example 5 shows, proving a similar statement for the parabolic case 
we must require some growth condition on the norms of products. 

Definition 3. We say that a gasket A : X m — > M n (K) is fast if there is a 
constant c > such that 

\\Au\\ > c\\Ar\\\\Aj\\. 
for every multi-indices I G X m and J G J m • X m . We call 

. \\A U \\ 

iei m \\At \\\\Aj\\ 
j&j m -x m 11 1111 11 

the coefficient of the gasket. 

Example 6. Consider the parabolic and hyperbolic gaskets of Example 5. 
Since 

\\A 1 (A 2 A k 1 \\ = ||^ +/£ ' +1 || = 1 + k + k' 

and inf fe fe /> 1 {(l + k + k')/{kk')} = 0, A is not fast. 

On the contrary, since any product of N = k + k! copies of B 12 is equal 
to B± and 

II rjk' r>k\\ H r>k+k' \\ rp \ Z7 Z7 II Tjk' II II 73k— 1 II 

||-t>i -t>i|| — H-Di || — r 2 k+2k' > t , 2k'r 2 k-2 — \\-E>i Ull-Di ||, 

then B is fast with c& > 1 . 

Example 7. The (parabolic) cubic gasket C 2 of Example 3 is fast. Indeed 
consider first J = 21L, with Cl = y ^ J , so that 



C, 



1 0\ fl l\ fa b\ _ f a + c b + d 
1 lHo 1 1 [c d ~ [a + 2c b + 2d 
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Clearly \\Cj\\ = max{a + 2c, b + 2d} < 2 max{a + c,b + d} and therefore 

\\MCj\\>±\\M\\\\Cj\\ 

for every M G SX 2 (N). TTie same argument applies to J = 12L. Since 
HCf'CiCaCf || = fc'(fc + 1), ||Cf'|| = fc' and HdCaCf || = 2A; + 1 it follows at 
once that in fact Cc 2 = 1/2. 

Example 8. A hyperbolic gasket is not necessarily fast. Consider for in- 
stance A:X 2 ^ GL 2 {Z[i]) with 

(2 0\ , f/3 
A ^={0 a )> A *={o 2 

where \a\,\(3\ = 1. A is hyperbolic since every Aj, \I\ = 2k, contains an 
entry with modulus 2 k and k! > k. On the other side A is not fast. Indeed 
WA^AlfW = 2- fc ||i4f||||Ai^|| for all k < k' and so 

iePjeJ* P/UPjII ~ °' 

Proposition 3. If A is a fast gasket with coefficient c, the gasket \A is a 
fast gasket for every A ^ with coefficient c/|A|. 

Proof. \\XAuW = \\W\AuW > clAH^m^H = ^A^H ||AA/|| □ 

Our main algebraic results on the exponent of a gasket are the following: 

Theorem 1 (Exponent of a fast gasket). If A : X m — > M n (K) is a fast gasket 
then < sa < oo. 

In the more detailed discussion in Section 3.2.2 we also show how to 
build explicitly two sequences of monotonically decreasing functions fA,k( s ) 
and gA,k( s ) such that, for all k from some k on, s^ G [^^(1), fa ^(1)] and 
I/4 ~ 9a]SS)\ — a ^°sk for some a > 0. These sequences will be used in 
Section 4 to evaluate analytical bounds for the exponents of a few semigroups. 

Theorem 2 (Alternate characterization of the exponent of a gasket). Let 
A : X m — > M n (K) be a gasket with s^ < 00. Then the function 



£a(s) = lim 

k— >oo 



En- 4 

\i\=k 



1 II 
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is a well defined log-convex positive strictly decreasing function and satisfies 
the following properties: 

1- U(P) = m; 

2- U(s) > 1 fors< s A ; 

3- UM = i; 

4- £a( s ) < 1 f or s > s A if A is a hyperbolic gasket; 

5- Ca( s ) — 1 f or s > s A if A is a fast parabolic gasket. 

Theorem 3 (Exponential growth of norms). Let A : X rn — > M n (K) be a 
gasket with s A < oo. Then 

r log N AM (k) 

hm = s A 

k^oo log K 

for every M G GL n {K). 

It is natural to ask whether this result can be made stronger, namely 
whether N A (k) x k SA . This leads to the following: 

Open Question 1. Based on numerical investigations, Boyd [Boy 82] con- 
jectured that, in case of the Apollonian gasket, the number N(k) of circles 
whose curvature is not larger than k (which, in our setting, corresponds to 
the number of Hirst matrices whose norm is not larger than k ) were only 
weakly asymptotic to k s , where s is the dimension of the Apollonian gasket, 
namely that 

N(k) ~ Ak s log* (k/B) 

for some A,B,t > 0. Recently Kontorovich and Oh [KOll] disproved this 
conjecture by showing that actually N{k) is strongly asymptotic to k s , namely 
N(k) x k s for all keN. 

We pose the following question: is N A (k) strongly asymptotic to k SA for 
every gasket A with a finite exponent? Or is there some extra condition that 
must be put on A to ensure this behaviour? 

Hausdorff dimension of limit sets of discrete subsemigroups of real 
and complex projective automorphisms. Recall that the projective 
space KP n ~ x is defined as the set of all 1-dimensional linear subspaces of 
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K n . We denote by [v] the 1-dimensional subspace of K containing v, so that 
[•] : K n — > KP n ~ l is a projection with kernel K \ {0}. Since linear maps 
preserve linear subspaces, every automorphism / : K n — > K n induces a 
projective automorphism ipf : PK n ~ l — > PK n ~ x defined by ^/([v]) = [/(f)] ■ 

We endow PK n ~ x with the round metric tensor, namely the metric with 
constant curvature equal to 1, and with the distance and measure induced 
by it. In every compact of every chart of KP n ~ x this distance is Lipschitz 
equivalent to the Euclidean distance in the chart and the measure is in the 
same measure class of the Lebesgue measure in the chart, so that we can 
forget about the definition and work directly with the Euclidean distance 
and Lebesgue measure in any chart. 

In this paper by iterated junction system (IFS) on a metric space (X, d) 
we mean a homomorphism F : X m — > C(X), where C(X) is the set of 

continuous maps from X into itself, such that the generators fi = F(i), i = 
1, ... ,m, are non-expansive maps, namely there exist Ki G (0, 1] such that 
d(fi( x ): fi(y)) < Kid(x, y) for all x,y G X. If Ki < 1 for all % then we say that 
the IFS is hyperbolic, otherwise that is 'parabolic. Recall that to a hyperbolic 
IFS F on X it corresponds a unique compact set Rp C X, which is also equal 
to the set of limit points of the orbit of almost every point x G X under the 
action of F, whose Hausdorff dimension, which we denote by dim H R F , is 
often non-integer and can be evaluated through some property of the /, e.g. 
their coefficients Ki or, if regular enough, their derivatives (e.g. see Chapter 
9 of [Fal90] for more details and examples). Finally, recall that the f\ satisfy 
the open set condition if there exists a relatively compact subset V C X such 
that V D \-i^Lifi(V). Note that this condition in particular implies that there 
is no relation between the fi. 

In a seminal paper D. Sullivan [Sul84] related the Hausdorff dimension 
of the limit set Rr C CP 1 of a (geometrically finite) Kleinian (i.e. discrete) 
subgroup T C PSL 2 (<C) to a critical exponent of T defined in the context 
of hyperbolic geometry (the recent results by Kontorovich and Oh quoted in 
Open Question 1 are in fact based on this result). In the same spirit here 
we relate (partly in form of conjecture) the exponent sa of a subsermgroup 
of SX^(R) or SL n (C) to the Hausdorff or Box dimensions of the limit set 
R A C KP n -\ K = R or C, of the subsemigroup of PSLf^R) or PSL n (C) 
naturally induced by A. This somehow complements recent results on the 
geometry of residual sets of real projective IFS by Barnsley and Vince [BV10] 
and on complex projective IFS by Vince [Vinl2]. 
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We denote by fi, . . . , f m linear volume-preserving automorphisms of K n , 
by Ai G SL^ (K) the matrices representing the fi in some coordinate system, 
by A the semigroup generated by the Ai and by ipi G PSL^ (K) the projective 
automorphism naturally induced by fi. 

Our main geometric results on the exponent of a gasket are the following: 

Theorem 4 (Semigroups of PSLf (R)). Let f u . . . , f m : R 2 ->■ M 2 6e suc/i 
i/ia£ £/ie induced maps tpi G PSXf (R) safe/?/ £/ie open sei condition with 
respect to some proper subset V C ILP 1 and £/ia£ ; m some affine chart if : 
M.P 1 — > R ; i/ie are contractions on (p(V) with respect to the Euclidean 
distance. Let Ra = (U\i\ = kipi(V)) be the corresponding residual set. 

Then 2 dim^ Ra = sa- 

The simplicity of the geometry of IRP 1 suggests that the theorem above 
can be extended to more general families of fi, leading to the following: 

Open Question 2. Let A C SLf(M) be a free finitely generated semigroup 
with finite exponent sa and ^/^ the corresponding free discrete subsemigroup 
of PSLf(M). Under which general assumptions is the Hausdorff dimension 
of the limit set of ^/a related to sa ? 

An identical theorem holds in the complex case: 

Theorem 5 (Semigroups of PSL 2 (C)). Let f 1: . . . , f m : C 2 ->■ C 2 be such 
that the induced maps ipi G PSL 2 (C) satisfy the open set condition with 
respect to some proper subset V C CP 1 and that, in some affine chart ip : 
CP 1 — ?■ R ; the ipi are contractions on (p(V) with respect to the Euclidean 
distance. Let Ra = H^ =1 (U\i\ = kipi(V)) be the corresponding residual set. 
Then 2 dim^ Ra = sa- 

Even in this case there is evidence that the claim of the theorem can be 
generalized further, at least to some parabolic case, namely when the ipi are 
not contractive but just non-expanding. This is, for example, the case for 
the Apollonian gasket as shown in the introduction. 

For n > 2 things get geometrically much more complicated for both R and 
C and we cannot claim any general result. In order to provide motivations 
for further studies of this case and as a source of examples of non-trivial fast 
gaskets, in Section 4 we introduce for free semigroups A C SX^(R) the notion 
of real and complex self -projective Sierpinski gaskets, the name being due to 
the fact that to such A's it correspond the same "cut-out" construction of 
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the standard Sierpinski gasket and of its multi-dimensional generalizations. 
In particular to each such gaskets it corresponds a curvilinear simplex Ta C 
KLpn-i invariant by the action induced on it by A and a compact invariant 
set Ra C Ta obtained by subtracting recursively curvilinear polyhedra from 

In particular in Section 4.3 we show that in PSL±(R), n > 3, there 
exist a smooth 1-parameter family of fast free semigroups A t , t G [l,oo), 
such that: 1) A t is a hyperbolic IFS for t G (1,4); 2) the extremal gasket 
A\ is the (parabolic, multi-dimensional generalization of the) cubic gasket; 
3) A 2 is the (multi-dimensional real self-pro jective generalization of the) 
standard Sierpinski gasket. Siilarly, in Section 4.2 we correspondingly show 
that in PSL2(C) there exist a smooth 1-parameter family of semigroups A t , 
t G [1/5, a], where a ~ 0.651, such that: 1) A t is a hyperbolic IFS for 
t G (1/5, a); 2) the extremal semigroup Ai/ 5 is the (parabolic) Apollonian 
gasket and Ai/ 2 is the (complex self-projective generalization of the) standard 
Sierpinski gasket. This shows that the cubic and Apollonian gaskets are both 
natural parabolic deformations of the standard Sierpinski gasket, the first in 
the real setting and the second in the complex setting. 

Finally, based on numerical and analytical results in Section 4.3, we claim 
the following about the box dimension of the residual set: 

Conjecture 1. Let A C SL^(M) be a real projective Sierpinski gasket and 
Ra its residual set. Then, under suitable natural assumptions, 

(n + 1) dims Ra > usa- 

3 Asymptotic growth of norms 

3.1 The case m = 1. 

Let us briefly discuss the exponent of a gasket in the trivial case m—1, since 
these results will be useful later in this section. 

Here X = N, A is generated by a single matrix A G M n (K), K = R or 
C, and A(k) = A k . In order for A to be a gasket then it is necessary and 
sufficient that A has an eigenvalue with modulus larger than 1 (in which case 
it will be a hyperbolic gasket) or that has a single eigenvalue of modulus 
1 and its spectral norm larger than 1 (in which case it will be a parabolic 
gasket). 
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If A is hyperbolic then \\A k \\ grows exponentially with k and so s = is 
the only exponent that can make the series (a( s ) — J2ken ll^ fc |l~ s divergent, 
i.e. sa = 0. 

If A is parabolic then its generator A has only eigenvalues of modulus 1 
and therefore its norm grows as some polynomial of degree d < n. Hence 
Ca( s ) diverges for s < 1/d and is finite for s > 1/d, i.e. sa = 1/d. 

Now consider the number A^(r) of powers of A whose norm is not 
larger than r. When \\A k \\ is a polynomial of order d their number grows 
as r 1 ^, so that the limit lim^oo log iV^r)/ log r = 1/d equals sa- When 
\\A k \\ grows exponentially then Na(t) grows logarithmically so that, again, 
lim^oo log A^4(r)/logr = equals s A - Incidentally, this proves Theorems 1 
and 3 for m = 1. 

3.2 The case m > 1. 

When there is more than one generator things are qualitatively different: 
the number of terms of order k (i.e. products of k generators) increases 
exponentially and there can be coexistence of polynomial and exponential 
growths of norms for terms of the same order. E.g. we already pointed out 
this behaviour for the C 2 semigroup in Example 7 and the same happens for 
all semigroups C n (see Section 4.3.2). 

Following Boyd's arguments in [Boy73b] we are able to find upper and 
lower bounds for Ca m ( s ) us i n g only the series of the norms of the "diagonal" 
terms A D M, D e V m , and of those "next-to-diagonal" ones AjM, J e J m . 
The key point of the next arguments is the following elementary recursive 
re-writing of the zeta function: 

CaM= iTTMp + Yl Ca Aj m{s). (2) 
Dev m " D " Jej m 

3.2.1 A fundamental set of inequalities 

Using (2) we can now write the following fundamental inequality: 
Proposition 4. Let A : l m — > M n (K) be a fast gasket with coefficient ca 
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and let J E J m ■ X m and I E I m . Then the following inequalities hold: 

n- s \\Aj\\- s Us) < ( Aaj (s) < c A s \\Aj\\- s Us) (3) 
Ca Ai (s) > va Ai (s) + n~ s n AAi (s)Ca(s) (4) 

CA Al (s) < VA Al (s) + C A S fl AA[ (s)Ca(s), (5) 

where va..(s) = tt— — — — and a a. As) = t— ; — — 

a m \ ) \\A D M \\ s ^ AmV ' ^ \\AjM\\ s 

Proof. The left side of (3) is a direct consequence of (1), its right side of the 
definition of fast gasket. The starting point to prove inequalities (4,5) is (2), 
from which we get 

Ca Ai (s) = v Aai (s) + Ca Aji (s)- 
Jej m 

Applying the right side of (3) to the summation above we get that 

E C>U,»< r s \\ Aj U s (a(s) = —VA Ai {s)^( S )- 

J&Jm j &Jm c A\\ A Jl\\ C A 

Hence (5) follows and analogously it is proven (4). □ 

Remark 1. The function v Am has the same complexity of the zeta functions 
of 1-gaskets. Indeed let Ai = (Ai) be the 1-gaskets generated by the m gener- 
ators of A. Then va m { s ) = Xa<i< m Ca( s )- ^ e ultimate idea of this section 
is to exploit this fact to find upper and lower bounds for ( A using the much 
simpler u A . 

Remark 2. For finite s, fi A and v A , converge or diverge together. Indeed if 
J = iD Ej m then J' = D E V m and 

11 4 " d |t < \\A W \\ < 7i||A||||^d||, 



n\\A 
so that 

(to — \)n~ s min || Ai\\^ s v A {s) < /i^(s) < (to — l)n s max ||A^ 1 || s z/^(s). 

l<i<m l<i<m 

The same is true for s — >• oo after some finite number of terms (those of "too 
small" norm) is removed from the two series. For similar reasons v Am and 
A* A M converge or diverge together with v A for every M E GL n {K). 
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Let us now examine closely the two inequalities (4,5) for 1 = 0. The first 
one becomes 

Ca(s) > »a(s) + n- s fi A (s)( A (s) (6) 

Since we are going to use this inequality to get lower bounds for ( A (s), we 
can proceed without loss of generality by assuming that ( A (s) < °°- Then 
we get that, for n~ s ji A (s) < 1, 



Us) > 



1-n s fi A (s)' 
Analogously the right one becomes 

Ca(s) < v A {s) + c A s fi A {s)( A {s). 



(7) 



(8) 



This time though, since we are going to use this inequality to provide upper 
bounds to the zeta function, we need to truncate the infinite series to a 
finite sum to ensure we are dealing with finite numbers. A natural recursive 
definition, inspired by the structure of (2), is 



Am k 



-°aJs) 



JeJ, 



e+i 



where v i is the necklace sum truncated at the order £. 
Proposition 5. Consider the sets V™ C X m recursively defined as 

'0 ~~ U Q i 

PI" = Wl) U v T+i-\j\ 



j 



e+i 



> 1. 



Then the following properties hold: 

1. Vf C V^ +1 for all £ > 0; 

2. U e > Vr = T m ; 

3- &>) = £76^11^11-. 
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Proof. 1. We prove this by induction assuming that V™ C V™ +1 for all 
k < I - 1. Now, if J e then either I eVf oi I = PJ with J G 
and P G In the first case I G P^ since V™ C P]^- In the second 

case we have that £ + 1 — | J\ < £ — 1 since every element of J7" m has at least 
rank two and therefore, by the inductive hypothesis, C 
Hence, by definition, 7 G P^+ 2 _|J| • J C P^+ r 

2. Notice first of all that every index / G X m either belongs to T> m (and 
so to some V™) or it can be factored out as a product I = J\ • ■ • with 
7j G JP™. This factorization simply consists in singling out the patterns 
of the form i\ ^ i 2 — 13 — • ■ • — i p inside I and is clearly unique. Let 
K = max!<j< fc I Ji\. By construction Jj G PJ£ for every 1 < % < k. Then 
'PTk+i ^ Ujej, m Vft-J contains all products PJj, 1 < i, j < k and similarly 
V™ K+1 contains all possible products of k of the Jj, namely I G V™ K+1 . 

3. Let us write Ca ( s ) = 21iegt\\ A i M \\~ 8 - Then if 1 G either 
||A/M|| _S appears in z/^ M (s), in which case I G P™ by definition, or in 

^A~a Af J '( s )' * n wn i cn case ^ = with K" G C^ +1 ~I J L This is exactly the rule 
that defines recursively the P™. Since we also have by definition of C^ M ( S ) 
that g° = V™ it follows that Q t = VJ l - □ 

Corollary 1. Let A be a fast gasket. Then the C.a m ( s ) satisfy the following 
properties: 

tAjs)<&)f°rall£>0-, (9) 

}im &„(s) = (10) 

n-\\Aj\\-'&{s) < Ca Aj (s) < cZ\\Aj\\-'(a(*)\ (11) 
< v Aai {s) + c-^ A/ (12) 

Proof. (9,10) are a direct consequence of points 1. and 2. of the previous 
proposition. (11) is a direct consequence of (1) (left) and of the definition of 
fast gasket (right). Using the rhs of (11) and then (9) we get that 

ci>) < E HA7M||-^ 1 -i J i( a )<ii Jf ( a )+ c -^(s)<sw 

from which (12) follows. □ 

Using the monotonicity in £ of v^ M (s) and A^ M (s) and setting M = t n 
we get that, in particular, 

&(*)<va(*) + ca 
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From this we deduce that, when /i^(s) < c s A , ( A (s) — \-c-^ A { s ) ^ or a ^ ^ anc ^ 
therefore, finally, 

In the next section we will use these bounds to build bounds for s A - 



3.2.2 The exponent of a fast gasket 

The idea of next theorem comes from the following observation, valid in the 
particular case when all generators A* of A have non-negative coefficients and 
norms larger than the maximum between 1 and the inverse of the coefficient 
C4 of the gasket. 

Consider (7). The function g^(s) = n~ s /j, A (s) is clearly monotonically 
decreasing with s, defined on (0, 00) and (7(0,00) = (0, 00). In particular 
the set 5u(s) < 1 is not empty and therefore (7) holds for s > s g , with 
s g = (?_ 4 1 (1). When we let s — > s+ the right hand side of (7) goes to infinity 
so that 0( s s) diverges, namely s A > s g . 

Now consider (13). The function /U(s) = c~ s fi A { s ) satisfies the same 
properties listed above for g A . Let Sf = /_^ 1 (1). Then /U(s) < 1 for s > Sf 
so that (7) holds. This proves that (a( s ) < 00 f° r & h s > s f, namely Sf > s^. 

This simple argument not only grants us that is finite but also provides 
for it non-trivial lower and upper bounds. In order to obtain a similar result 
in full generality we need to refine (13). This will lead us to refine also (7) 
and to generate a pair of sequences converging to s A from the left and from 
the right at logarithmic speed. The idea is to apply over and over recursively 
first the inequalities (4,5) and then the inequality (3) to the truncated zeta 
function. 

The starting point is the sequence of sets of multi-indices Q^,k built as 
follows. We define Q\ k = J m ■ Then we consider the sets recursively defined 

as 

Q% k = {J e Q%1 1 > k) \Jj m ■ {J e J m \ J e Q%1 \\Aj\\ < k} 
with t>\. 

Proposition 6. For every gasket A and every k > there exists a I such 
thatQ% k = Q%l 
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Proof. The only thing that the recursive algorithm does is replacing the 
indices of Q e A k corresponding to matrices with norm not larger than k with 
indices of higher order. The set Q^l thus obtained can still contain indices 
corresponding to matrices with norm not larger than k but in a finite number 
of steps all such indices will disappear because, by definition, there is only a 
finite amount of them. Hence there exists a finite £ such that \\Aj\\ > k for 
all / G Q e A ,k- Then the algorithm leaves Q A ,k unchanged. □ 

Definition 4. We use the notation Q.A,k = Q-A,k an d> correspondingly, 

fU s ) = c A E ii^ii-, fAjk{s) = cx J2 ( 14 ) 

9A,k(s)=n- s J2 ll^ll"*' 9Ajk{s)=*T' J2 \\M\~'- (15) 

Note that fA,o( s ) — an d 9A,o( s ) = 9a( s )- Next proposition though 

shows that, for k big enough, fA,k and gA,k have a nicer behaviour than /a 
and gA do for a general A. 

Proposition 7. For every fast gasket A there exists a k and a ^A,k < oo such 
that fA,k an d gA,k are strictly decreasing continuous functions of s defined in 
(lA,k, oo) and with image (0, oo) for all k > k. 

Proof. By construction fA,k(s) and gA,k( s ) are proportional to each other 
and, respectively with constants c~^ and n~ s , to the sum of a finite number 
of functions jjL^ A defined as the series [Ia Ai fr° m which all terms with norm 
smaller than k have been subtracted. 
By Remarks 1 and 2 then 

n- s min HAH -5 ^ (s) < 1$ (s) < (m - l)n s max II A" 1 !! 5 V zA fe) (s), 

KKm " 11 Aa i K ' ~ r>AA / v ' ~ v ' Ki<m" Aa i 

leg 

where Q C X m is some finite set of indices and the series i/^ is equal to va Ai 
minus those terms with indices D G T> rn such that ||^4j/|| < k for all J G J rn 
with J' = D. In particular then also fA,k and gA,k(s) are bounded by a finite 
sum of series (s) and so, by the case m — 1 of Theorem 1 discussed in 
Section 3.1, they are finite on some connected non-empty interval (7.4^, 00). 
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Now let k = n 2 max!<j< m H^ 1 )! max i<i<m ||A||- F° r every k > k we have 
that, if the index DI appears in the series Vj^ and J is one of the indices 
such that \\Aji\\ > k, 



namely 



n max \\A i 1 || n 2 max \\A i 1 || max 

l<i<m , KKm l<i<m _ 
< — — < 1. 



\\Adi\\ ~ k 
This means that all summands in f Ak and g A ^ k are s-powers of numbers 
uniformly bounded from above by some number strictly smaller than 1 and 
therefore f Ajk and g Atk are strictly decreasing with s and their image equals 
(0,oo). □ 

Theorem 1. Let A be a fast gasket of m matrices. Then < s A < oo and 
both s g> k = 9a\{1) an d s f,k = /.Afe(l) are uniquely defined and converge, re- 
spectively from left and right, to s A as k — > oo with speed at least logarithmic. 

Proof. We start by showing that 

(a(s) > h Ajk (s) + g A ,k(s)( A (s) (16) 

and 

Ca(s) < h A , k (s) + f A , k (s)( A (s) (17) 

for every k, where h Ak is some positive continuous function that plays no 
role here. We prove it in detail for the most complicated case, namely the 
second one, when we need to pass through the partial sums ( A . First we 
notice that (12) writes as 

C AA ^)<h AAi , k + f AAi , k (s)C A (s) (18) 

after putting H Aa fc (s) = v Aai (s). Now we write the definition of ( Aa (s) 
splitting the last term in two 

<S*w = <w + E <£; w w + E <5£ w m a*) 

||Aj||<fc HA/||>fc 
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and apply (12) to the first term in ( and (18) to the second, getting 

\\Aj\\<k \\Aj\\>k 
= h\ Ai , k {s)+C A S \\M-'&(*)=h\ AiJc (8)+fX AlJt (8)&(8), 

where we set h\ k (s) = h° A Js) + E/e^» va Aji {s). 

1 1 \\Aj\\<k JI 

By repeating recursively this procedure we get that 

< h\ k {s) + f r Ak (s)Us) 
for all r > 0. Since the Q\ k stabilize into the Q A ,k, we proved that 

( A (s) < h Atk {s) + f A , k (s)( A (s). 

for every £ and therefore (17) follows; (16) is proved analogously. 

By Proposition 7, for k big enough the points Sf >k = f A \(l) and s 9jk = 
g A \(l) are uniquely defined and strictly between and oo. An argument 
analog to the one at the beginning of the section shows that these inequalities 
imply that s 9jk < s A < Sf jk for every such k. In particular < s A < oo since 
every Sf tk is finite and every s 9tk is larger than 0. 

To prove the last part of the theorem we point out that g A , k (s) = c s A n~ s f A , k (s) 
and that, for any s' > 0, 

/*,*(«) = ^ E \\ A A\- S > (c A ky ,n ^ J2 \\ A A- s - s ' = (c A ky'f Atk (s>). 



Hence 



so that 



1 = 9AA s 9,k) = c S9 - k n s °' k f A>k (s g>k ) > 
> c s X' k n- s ^{c A k) s ^- s ^f Ak {s Lk ) = c'jfn-'-xk 8 **- 9 "* 



> s/,fc(log k + log c A ) - s g>k (log k + log n) 

and finally 

/log A; + log n \ logn-logc^ 

< s Lk - s g , k < s g , k — — — - 1 < s A 



log k + log C4 / " log k + log c A 

□ 
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3.2.3 An alternate characterization of s A 



The exponent s A of a m-gasket A can also be extracted from the asymptotics 
of the partial sums of the ||Ar||~ s over same-rank multi-indices, namely from 
the sequence of functions (A,k( s ) = Yliex m II^MI _S - 

Lemma 1. The sequence of analytical log-convex monotonically decreasing 
functions C](k( s ) converges pointwise for every s G [0, oo) to a bounded con- 
tinuous log-convex monotonically decreasing function £a( s )- 

Proof. The (^l(s) are analytical log-convex monotonically decreasing func- 
tions because every summand ||A/||~ S of the satisfies those properties 
and so does every finite or infinite (converging) sum of them and positive 
power. 

In order to prove the convergence of the sequence we can replace with- 
out loss of generality in the expression of the (by abuse of notation we 
will denote the new functions still by C^fe) the max norm with any sub- 



multiplicative norm 



and notice that, since II A 



u\ 



then 



(A,k+k'(s) > CA,k(s)CA,k'(s)- From this follows that, for every s, the sequence 
Cj(k( s ) can have only one accumulation point and this point must be equal to 



A/ k ( 
su PfcGN SA,k\ S 

other elements of the sequence are not smaller than (s 



The main point is that for every element (J1°,(s) almost all 



i/fc 



0. Indeed if k = Nk then immediately C/ui( s ) > [(Ak ( s 



l/feo 
A,k 

e for every e > 

l/k _ ^l/ko 



A,k \ s )i 



while if k 

that 



Nk + £, with 1 < t < k - 1, then ( A:k (s) > (% ko (s)( A A s )i so 



£&(*)> cZK'W&t 

Since there are only a finite number of possible values of 
we can find a N big enough s.t. both 



Nk +l 



c 



1/fco, 
A,k { 



l+t/(Nk ) ujvfe 0+ ^/ ^ 

SA.e V s ) 



A,e 

for every e' > 



art*) 



A/ko , 

^A,k y 

< e' hold for all t. Hence 



l+£/(JVfc ) 



^A,k 



s) 



< e' and 



& > - e> (C%»{ 8 ) + l-e')> Cfklis) 



for small enough e' . 



That £a( s ) = su PfceN C A k( s ) be finite for all s is clear from the fact that 



all Ci(fc( s ) are positive decreasing functions bounded by CXki^) = m - 
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.i/fc. 



□ 



Theorem 2. The function £4 satisfies the following properties: 
1. £ A (s) > 1 for s < s A ; 

2- UM = i; 

3. £a( s ) < 1 f or s > s A if A is a hyperbolic gasket; 
4- £,a( s ) — 1 f or s > sa if A is a parabolic fast gasket. 

Proof. Directly from the n-th root test we get that £a( s ) > 1 for s < s.4 and 
£a{ s ) < 1 f° r s > S A> so that in particular £a( s a) = 1- 

Assume first that A is hyperbolic, so that there exist constants a > 1 
and K > s.t. ||A 7 || > ifaM for every I e X m . Hence 

namely for every e > we can find a o/ > 1, with \a — a'\ < e, and a 
k > such that (In C fe 1/fe (s))' < - In a' for all k > k. Since log a Cj^M = 
this means that, for every k G N and s > 0, ln(^/ fc (s_4 + s) < —sin a' and 
ln(^/ fc (s.4 — s) > sin a', namely Cfc^( s ^) > ( a/ ) s > 1 at the left of Sa and 
Cy fc (s^) < (a')~ s < 1 at its right. Since this is true for almost all k, the 
same properties hold for £4. 

Assume now that A is fast parabolic and that sa < 00 (e.g. in case that 
A is fast). In this case the sequence = min|/| =fc {||A/||} grows polynomially 
and therefore, for s > sa, 

1 > Cl /k (s) > a- k s/k = (a-JY S ^ 1, i.e. Us) = }™ (1%) = 1. 

Proving that £a(s) > 1 for s < sa requires much more work. We start by 
observing that, analogously to (3) and (2), we have the inequality 

C^M > ^fA^^M (20) 
and we can re-write CA Al ,k as follows: 

Cu jlfc ( s ) = Yl Ca Adi ,o( s )+ Yl (A AjI ,k-\J\(s). (21) 
dgx>™ Jej™ 

\D\=k 
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Applying (20) to (21) we get that 

Cu 7 »> E C^ 7 ,o(s)+n- s Uji\\-{A,k-\J\(s). 
Dev™ Jej™ 

\D\=k 

We set w°f (s) = YsJ^J? n ~\\ A Ji\\~ s , so that 

\J\=j 

Cu,»> E w A Al J s )(A, k -j(s). (22) 

2<j<fc 

Now we split the second summation in (21) as we did in (19) and obtain 

(A Al , k (s)= (Aa di ,o( S )+ E (A AjI ,k-\J\(s) + E (A AjI ,k-\J\(s). 

Dev™ JeJ^ JeJ£> 

\D\=k \\Aj\\<* \\Aj\\>* 

Now we apply recursively (22) to the first term in the splitting and (20) to 
the second, obtaining 

(A Al , k (s)> E (A AjI ,k-\j\(s) + (A AjI ,k-\j\(s) > 

JeJ™ Jejjj 1 

\\Aj\\<K \\Aj\\>k 



> E E ™A AjI Js)CA, k -A*)+n- s E \\Aji\\- s Ca^ui(s) = 
2<j<k Jej™ Jej™ 

\\Aj\\<K \\Aj\\>K 

= E w A Al A s )Uk-j(s), 

2<i<fe 

where _ 

«il» = E ^ AjI ,M+n- J2 UA-- 

\J\=j \J\=3 
\\Aj\\< K \\Aj\\>k 

Since there are only finitely many terms in A with norm not larger than k, by 
repeating recursively this procedure, the w l j£ A K stabilize to some functions 

w a Ai ,k such that Ca Ai A s ) > E 2 <i<fc^ AjlK ( s )CAfc-i( s )- For / = we get 

(Ak( s ) ^ E w A,MCA,k-j{s). 
2<i<fc 
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Now consider the polynomials pfc(x) = x k — J2 2 <j<k W A S s ) xk _J - By Descartes' 
rule of signs they all have a single positive root. Moreover for k big enough 
this root is larger than 1. Indeed by comparing the recursive algorithm that 
generate the w A K with the one generating g^ K it is clear that the function 

W A Al A s ) = T,2<j<k w A Al A s "> is e( l ual to the truncation of the series g AAj , K 
to the matrices Ai with |/| < k. In particular this means that when s < s A 
there is some k large enough that W\ K (s) > 1. Then pi(l) < and 
therefore the positive root a of pj. is larger than 1. Hence 

p = min {{ A ,j(s)<T~ 3 } = inf {CAj( s ) a ~ 3 h 

0<j<k 0<?<oc 

which follows by induction as a consequence of the following observation: 
(ami(s) > E w a Ai , k (s)(aMi-M > ^ E < Al ^ s V +1 ~ j > ^ S+1 - 

2<j<fe 2<j"<fc 

Hence lim^oo C.Y j( s ) > o" > 1 for every s > s A . □ 
Example 9. Let Mi, . . . , M m be upper triangular matrices of the form 



Mi 

and assume that 



1 



max < -. 

<i<m 1 - maxi<i< m {l - ai) 



It is easy to prove by induction that under this assumption the non-zero off- 
diagonal term never gets larger than 1, so that ||Mj|| = 1 for every I G X m . 

Now consider the gasket A generated by Ai = PiM i} > 1. By the 
observation above for every I — i ± . . . i k we have that \\Aj\\ — p^ • • ■ p ik and 
therefore 

Cu = E \\ A *\\- = E • ■ • pi s = 0°r s + ---+p~ m s ) k - 



lei™ igi™ 



Since A is clearly a hyperbolic gasket, by Theorem 2 its exponent s A is the 
unique solution of the equation p± s + • • • + p~ s = 1. Similar but more com- 
plicated conditions can be found for upper triangular matrices in higher di- 
mension (see Section 4-3.1 for a similar case with 3x3 matrices). 
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3.2.4 Norm asymptotics of fast gaskets 

We can use the previous section's results to prove Theorem 2. 

Lemma 2. Let A : X m — > M n (K) be a semigroup homomorphism and M G 
GL n (K). Then 

N ^ > N - (^m) ' < 23 » 

Proof. Since ||Afcf|| < n||A|| ||M||, we have that < ^ =^ < 
r,naindy{A J |p J ||<^}c{^jM|p J M||<r} " □ 

Theorem 3. Let A be a hyperbolic or fast parabolic gasket. Then 

v log^ M (r) 

urn = sa. 

r— >oo log r 

for every M G GL n (K). 

Proof. Since ||^4||/(n|| _1 1|) < ||AM|| < n||A||||M|| we can prove the theo- 
rem without loss of generality for M = t n . 
log N A (r) 

limsup < Let s > s^. Then 

r^oo log r 



>Ca(s)> Yl WM\~ S > E r- s = N A (r)r- s , 

so that 



oo 

\Ai\\<r \\Ar\\<r 



s | logq(g) > iggjW 

log r log r 

and therefore limsup r ^ 00 log 1 ^, ( - r - ) < s. Since this is true for every s > s A it 
follows at once that limsup,,^^ log 1 ^ r ^ < s^. 
log N A (r) 

lim inf > s A . From the elementary observation that 

r— >oo log V 

{A:\WAjW <r,IeI m }D |J {AuIWAjjW <r,Iel m } 
and using (23) we get that, for every k G N, 
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Assume now s > s^. Since A is, by hypothesis, either a hyperbolic or a 
fast parabolic gasket, by the definition of gasket and Theorem 2 we can always 



choose a ko such that ||Ar|| > l/ n f° r l-^l — and 
Now set a m 

s.t. N^(r ) > and set r\ = a^To and rj = a, 
proof of Theorem 2 we have that, by induction, 



> n s 



fimin/pT™ IIAHI and au = nmax lG r^ \\At\\, let r > be 

- fc " feo " 

' -1 ri, % > 2. Similarly to the 



M 



def 



min N^(ro)r 

re[ro,ri] 



inf A^4 (r )r~ 

rG[r ,oo] 



Indeed note first of all that lim^^ = oo since we chose &o so that a m > 1. 
Assume now that N^(r) > Mr s in [ro,Tj] and let r e [r,, rj+i]. Then for 
every J e X£J we have that r/(n||A j| 



> 



G [ro, rjj since 



n\\A. t 



> 



and therefore 



n||A„ 



> 



n||Aj 



> Mr s . 



Hence it follows at once that 

logiVA(r) 
logr 



logM 

> -j + s 

logr 



and therefore liminf^^^ l ° 9 , N ^ r ^ > s. Since this is true for all s < it 



follows that liminL 



■ r^oo logr 



> SA- 



□ 



4 Hausdorff dimension of limit sets of dis- 
crete subsemigroups of real and complex 
projective automorphisms. 

In this section we show how the exponent of a free finitely generated semi- 
group A C SL n (K) (resp. A C SL±(K) if n = 2n' and if = R) is 
sometimes related to the Hausdorff dimension of the set of limit points of 
a generic orbit in KP n ~ x of the subsemigroup ^(A) C PSL n (K) (resp. 



32 



ty(A) C PSL^ n ,(R)) naturally induced by A (equivalently, to the residual 
set of the IFS corresponding to ty(A)). 

Throughout this section we will make analytical and numerical evalu- 
ations of the exponent of several semigroups. The analytical bounds are 
obtained via the functions /iA,k defined in Section 3.2.2 thanks to Theorem 
1. The numerical ones are obtained by evaluating numerically the function 
NA{k) and interpolating the curve log NA(k) with respect to log A; thanks to 
Theorem 3. Since A is a gasket and NA{k) is integer- valued, a computer 
program can evaluate exactly its values, the only constraint coming from the 
running time which increases exponentially with k. All calculations were 
done with Xeon 2MHz CPUs under Linux. 

4.1 n = 2,K = R 

Let {/i, • • • , / m } De a free set of linear automorphisms of R 2 preserving a 
volume 2-form modulo sign. With respect to any frame £ = {ei,e 2 }, these 
automorphisms are represented by matrices Ai G SXf (IR). We denote by A 
the semigroup generated by the Ai and by ipj G PSLf(M) the automorphism 
of 1LP 1 ~ S 1 naturally induced by //, I G X 2 . The similarity between the 
characterization of the exponent sa given in Theorem 2 and the formula for 
the Hausdorff dimension of a 1-dimensional IFS given in [Fal90] (Theorem 9.9, 
p. 126) suggests the following claim: 

Theorem 4. Assume that the fa have all real distinct eigenvalues and that 
there exists some proper open set V C M.P 1 invariant by the ipi and such 
that, for some affine chart (p : 1LP 1 — > R, the ipi are contractions on <f(V) 
with respect to the Euclidean distance and satisfy the "open set condition" 
4>i{V) n i> 2 (V) = 0. Let R A = C\^ =1 (U^^/iy)) be the corresponding 
residual set. Then 2 dim H Ra = sa ■ 

Proof. Let (x, y) and [x : y] the affine and homogeneous coordinates asso- 
ciated to <f and assume, for the argument's sake, that ip([x : y]) = x/y. 
Under the hypotheses every ipj has two fixed points, an attractive one aj 
and a repulsive one rj. Since we are assuming V to be invariant under the 

ipi, then it must happen that a/ G V for all / G X m . Let A/ = ^ be 

the matrix representing the // in the affine coordinates relative to if. Let 
= (^o^o^^x) be the coordinate expression of ipi in the chart ip. Then 
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a direct calculation shows that 



aip + (3 



det Aj 



(recall that det Aj = ±1). Now let <f(V) = [ipi, (p 2 ] and set ip Ti 
Then ||A/||/2 < fnp m + v < \AA\ and therefore 



(y?i + y? 2 )/2. 



1 



Hence the limit 



lim 

k— >oo 



< 



< 



\\M\ r 



\I\=k 



converges iff it converges the limit 



lim 

k— ¥oo 



\i\=k 



By the result on 1-dimensional contractions quoted above, the exponent sep- 
arating the values of s for which the first limit diverge from those for which 
it converges is exactly dim# Ra- By Theorem 2 this means exactly that 
2 &iva H Ra = sa- □ 

An interesting consequence of the previous theorem is the following con- 
straint posed by geometry to the (algebraic) exponent of the semigroups 
satisfying its conditions: 



Corollary 2. Let A C SLf' 



j 2 v^; be a semigroup satisfying the conditions of 
the theorem above. Then sa < 2. 



Proof. This is a direct consequence of the trivial fact that the Hausdorff 
dimension of a subset of R cannot be bigger than 1. □ 



4.1.1 Matrices with non-negative entries 

The semigroup S'L^ : (]R + ) is a source of several interesting semigroups that 
satisfy the hypothesis of Theorem 4. Indeed every linear automorphism 
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/ = A % -ti ® preserves the positive cone C{£) over £ and therefore the in- 
duced projective automorphism ipf preserves the closed segment S(£) C RP 1 . 

In this simple setting there is an easy sufficient condition to determine 
whether a gasket is fast: 

Proposition 8. Let A : I m — > SLf(EL + ) be a gasket and suppose that all 
products Aj, \I\ = 2, have no entry equal to zero. Then A is fast. 

f D Q\ I Qj b\ 

Proof. Let A 12 = ( J and A K — ( ^ j , so that 



A\2K — 



pa + qc pb + qd 
ra + sc rb + bd 



Assume, for the argument's sake, that ||Ai2A:|| = ra + sb. Then, since 

minjjo, q] 

pa + qc> -{ra + so), 

maxjr, s} 

we have that, for every M G M 2 (1R + ), 

\\MA 12K \\ > \\M\\{pa + qc) > mm{ /' Q \ \\M\\ \\A 12K \\. 

maxjr, s} 

By repeating this argument for every index of order 2 and denoting by c 
the smallest of these coefficients, we have that ||MAi2A-|| > c||M|| || ^4i2it || for 
every M G M 2 (IR + ). In particular then A is a fast gasket with coefficient not 
smaller than c. □ 

Example 10. Let £ = {e±, e 2 } be a frame on IR 2 and / 12 defined by 

/i(ei) = ei + e 2 , /i(e 2 ) = e 2 ; /2(ei) = 2e x + e 2 , / 2 (e 2 ) = e 2 . 

With respect to £ the f\ are represented by the matrices 



Pi 



i or 2 u o 



The semigroup F = (F 1 ,F 2 ) C SL~(N) is free because if F r G A, I G X 2 ; 
t/ien i/ie entries in Fj 's lower row are equal to the entries in the upper row 
of the matrix Fp and according to whether the upper left entry of Fj is larger 
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where fk, k > 1, is the Fibonacci sequence 0, 1, 1, 2, 3, 5, . . . 



or smaller than its lower left entry we get whether I = 21' or I = II'. 
Proceeding recursively this way we see that there is no other index J ^ I s.t. 
Fj — Fj. In particular then F is a gasket. Moreover F is hyperbolic: indeed 
clearly > \\Fl\ since F 2 has no entry smaller than the corresponding 
entry of F 1; and ~ g k , where g is the golden ratio, because clearly F* = 

fk+2 fk+1 
fk+1 fk 

Finally, F is fast (with coefficient not smaller than 1/3) by the previous 
proposition. 

A direct calculation shows that i[) 1>2 are not contractions over S(£) but 
they are so over the smaller set S{£'), with 

£' = {e\ = (1 + V3, 2), e' 2 = (1 + Vs, 1)}. 

Let [x : y] be homogeneous coordinates on KLP 1 corresponding to £' . In the 
canonical chart ip = x/y the maps ipi induced by f\ write as 

i ( \ V + l i ( \ 2 V + 1 



if ip 

which reveals that this example coincides with Example 9. 8 of [FaWO], coming 
from the theory of continued fractions. 

To obtain analytical bounds for s F we can use Theorem 2. Since both 
generators have an eigenvalue larger than 1 the norms of the terms F\F 2 
and F 2 F± grow exponentially, so that we can get a good approximation of 
^F,e by truncating the sums after just a few terms. By considering only 
the terms with k < 10 in fi F and solving the equation Hf,o( s ) — 2 s in 
this approximation we get sp > -51, with a relative error of about 6% on 
the more precise estimate sp > .54 obtained by considering k > 20. Since 
c = 1/3, the first fiF,e we can get upper bounds is fiF,3- Here we just mention 
that from /if,8, considering the first 30 summands of all series that appear 
in its expression, we get 0.95 < s F < 1.76. In terms of the dimension of 
R F this translates in 0.474 < dim H R F < 0.877. By evaluating N F (k) for 
k = 2 P , 1 < p < 28, (taking about 20 minutes of CPU time) we get the 
estimate sf — 1.062 (see Table 3 for the corresponding values of Np), with 
a (heuristic) error of 2 on the last digit. This corresponds to the well-known 
fact dim H R F ^ 0.531 . 

Example 11. Consider now 

/i(ei) = ei, /i(e 2 ) = e ± + e 2 ; /2(ei) = e x + e 2 , f 2 (e 2 ) = e 2 . 
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The corresponding matrices ( with respect to S) 



Ci- (J lj,c 2 - (J \ 

generate the semigroup C 2 C SX 2 (N) we already met in Examples 3 and 7. 
In particular we already know that C2 is a parabolic fast gasket with coeffi- 
cient c = 1/2. A direct check shows that the slowest and fastest growths with 
respect to the order k of the multi-index I of Ci € A correspond respectively 
to the pure powers C\ , for which \\C\\ = k, and to the "cyclic" products 
CiC i+1 C i+ 2- ■ -Ci+k-i, for which ||CiC i+ iC i+2 • • • C i+fc _i|| ~ g k , where the 
sums in the indices are meant "modulo 2" in the sense that 3 means 1, 4 
means 2 and so on. The reason why the golden ration g appears is that, 
similarly to the previous case, 



n n n n ( f k + 2 /w I 

O1O2O1 ■ ■ ■ — 1 , , 1 

J h J 



for k odd while if k is even the two rows get exchanged and analogously for 
the cyclic products beginning by C 2 - 

In the affine chart ip : [x : y] — >■ x/(x + y) the maps ipi induced by the f\ 
write as 

1 + (p 1 — <p 

and the segment S(S) maps into [0, 1] . Note that this choice of chart corre- 
sponds to writing e\ = e[ and e 2 = e[ + e' 2 , expressing the f\ with respect to 
£' = {e'^e'2} and using the canonical chart y' — 1 for the corresponding ho- 
mogeneous coordinates [x' : y'\. In terms of the semigroup, this corresponds 

to the adjunction via the matrix M = \ \). A direct calculation shows 



that the |^(y?)| < 1 on [0, 1], with the equal sign holding at for tpi and at 
1 for ip2, namely the IFS {1^1,1^2} is parabolic. 

Evaluating the Hausdorff dimension of the limit set Rc 2 of a point w G 
(0, 1) under the action of C 2 is nevertheless an easy task. Indeed, since 
^((0,1)) = (0,1/2) andi/j 2 ((0,l)) = (V 2 ,!), the fa, |/| = k, subdivide (0,1) 
into 2 k disjoint segments di = ipi(0,l) in such a way that U\i\=kdi = [0,1]. 
Moreover the length of these segments goes to zero for k — > 00. Indeed if 
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Cl = (c d) e C2 tkm ^ /(V?) = (a+tc%lc+d and therefore 



c + d a + b 



(a + b)(c + dY 



from which we get 

< \di\ < 



4||C/|| 2 - 1 J| - \\d\\ 

Hence the orbit under the ipi of every w G (0, 1) is dense in (0, 1) and there- 
fore dim H R c . 2 = 1 . Note that, since this IFS is not hyperbolic, Theorem 4 
does not apply to it. 

Let us get analytical bounds for via Hc 2 ,e- Unlike the previous example, 
the presence of parabolic elements in the semigroup does not allow to truncate 
the series to just a few terms because of its very slow convergence. Since 
\\CiCi +1 \\ — k + 1 we get easily that 

oo oo 

^c 2 ,o(s) = £ UCjII- = 2£ I|CiC*||- = 2£V* = 2(C(*) - 1), 
jej 2 k=i k=2 

where £(s) is the Riemann's zeta function. The solution of /ic 2 ,o( s ) — 2 s 
gives us the bound sc 2 > 1.54. The first upper bound can be gotten from 
k L c 2 ,2, obtained by replacing the two terms of norm 2 in /ic 2 ,o( s ) , namely 
1 1 C12 1| s and ||C2i||~ S ; with, respectively, fic 2 A 12 ( s ) an d fic 2 A 21 (s). A direct 
calculation shows that 

oo oo _ 

Vc 2 a 12 (s) = £(2A; + 1)- + E k ~ S = 2 ~^(s, + " 1 " 2 " " 3 ^> 

k=2 k=4 

where C(-M) the Hurwitz zeta function, and by symmetry we know that 
^c 2 Ai2 = ^c 2 a 21 ■ Hence 

HcM = 2 ( 2 C(*) + 2~ s C(s, \) ~ 2 - 2 1 - - 3-^ 

which gives the bounds 1.7 < sc 2 < 3.93 as solutions of nc 2 ,2( s ) = 2 s and 
l^c 2 ,2(s) = 2~ s . yl numerical evaluation of Nc 2 {k) for k = 2 P , 1 < p < 17, 
(taking about 1 hour of CPU time, see Table 3 for the evaluated values) gives 
Sc 2 = 2.0001 with a (heuristic) error of 1 on the last digit. This and the 
evaluation of dim h Rc 2 above strongly suggest that sc 2 = 2. 
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It is interesting to consider the following generalization of the previous 
example, namely the free semigroups C 2 , a C SL 2 (N) generated by 




a 
l/a l/a 



) 




l/a l/a 
a 



In this case, in the same framework used above, 



a 2 + (2 - a 2 )v? 



1 + (a 2 - l)y? 

2 + (a 2 -2)^' 



A direct check shows that, for every fixed a G (1, 2), the ip^ a are contractions 
on the invariant interval [0, 1] and that they satisfy the open set condition 
with respect to it. Let Rc 2a be the limit set of the orbit of any point 



4 we get the following: 

Proposition 9. s c . 2>a = 2 for every a G (1,2). 

Remark 3. The restriction on the possible values of a looks more like an 
artificial effect of a poor choice for the distance function rather than a true 
property of the semigroups. We believe that by choosing a ad-hoc metric and 
maybe slightly modifying the argument the proposition above can be extended 
to the half-line [l,oo). 



Now we consider the case of 2 x 2 complex matrices. The corresponding 
projective space is the Riemann sphere CP 1 , namely the complex plane with 
the addition of a point at infinity. 

The geometry of Kleinian groups, namely of discrete subgroups of the 
Mobius group PSX 2 (C), is known to be extremely rich and is presumably 
even richer in case of Kleinian subsemzgroups. A study of such semigroups in 
a general setting is therefore way beyond the scope of the present paper. Here 
we rather state first a somehow general property analogous to the real case 
above and then, as a source of examples, focus our attention on a particular 
but interesting case that we refer to as complex projective Sierpinski gaskets. 




4.2 n = 2,K = £ 
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Theorem 5. Let /j be m linear automorphisms of C 2 represented in coor- 
dinates by the matrices Ai e SL 2 (C) and let ipi be the induced elements in 
PSL 2 (C). Assume that there exists some proper open set V C CP 1 invari- 
ant by the tpi and such that, for some affine chart <p : CP 1 — > C, the ipi are 
contractions on <p(V) with respect to the Euclidean distance and satisfy the 
"open set condition" ^i(V) H ip 2 (V) = 0. Let R A = n^ =1 (\J\i\= k i/>i(V)) be 
the corresponding residual set. Then 2dimnRA = sa- 

The proof of this theorem is almost verbatim the same of Theorem 4 and, 
correspondingly, we have the following corollary: 

Corollary 3. Under the hypotheses of the previous theorem, 2 < sa < 4. 

Definition 5. Let /i,/2,/3 be volume-preserving linear automorphisms of 
C 2 with real spectrum and denote by A ± ,A 2 ,A 3 e SL 2 (C) the corresponding 
matrices with respect to some coordinate system and by ip 2 , i>3 £ PSL 2 (C) 
their corresponding projective automorphisms acting on the Riemann sphere. 
Let [ei] E CP 1 be a fixed point for ipi corresponding to the largest eigen- 
value A > 1 o//j. We say that the semigroup F generated by the fi (or, 
equivalently, the semigroup A generated by the A; L ) is a complex projective 
Sierpinski gasket if the following conditions are satisfied: 

L [fi( e i)] = [fj( e 'd] f° r ever v hj with i ^ j; 

2. the circle Tk passing through [e*], [e^] (where i,j,k is a permutation of 
1,2,3) and [fi(ej)] is invariant under both fi and fj, i ^ j; 

3- and [fki&j)\ belong to the same connected component of€,P l \Yk 
for every permutation k) of (1, 2, 3). 

4- the circles T k do not intersect in the interior of the curvilinear triangle 
Ta C CP 1 having as vertices the [e*] and as sides the segments of the 
Tk with vertices the points [ej] and [ej] containing the point [f%(ej)]. 

By construction every such gasket A is free and satisfies the open set 
condition with respect to the interior of Ta- Since the Mobius group PSL 2 (C) 
is transitive on triples of distinct points, we assume without loss of generality 
in the rest of this section that Ta has vertices [ei] = [1:1], [e 2 ] — [i : 1], 
[e 3 ] = [—1:1] with respect to homogeneous coordinates [z : w] and use the 
affine chart w — 1 with complex coordinate z = x + iy for all calculations. 
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Proposition 10. Let /i,/2,/3 be volume-preserving linear automorphisms 
with real spectrum having respectively e\ = (1, 1), C2 = (i, 1), e 3 = (—1, 1) as 
eigenvectors corresponding to the largest eigenvalue and assume that 

^i(N) = V^Oi]) =u + iv, 
^(N) = ^i([e 2 ]) = is, 
.^(N) = ^2([e 3 ]) = -w + iv. 

A necessary condition for fi, / 2 , f 3 to generate a Sierpinski gasket symmetric 
with respect to the imaginary axes, namely such that fi(z) = f2(—z) an d 
f 3 (z) = f 3 {—z), is that ipii\e 3 ]) G V , where Y is the circle 

x 2 + y 2 -x(l- s 2 ) - s 2 = 0. (24) 

For s = the condition is sufficient for u e [1/5, a], where a ~ 0.651. 

Proof. A long but straightforward direct calculation shows that condition 
(24) is the only one coming from imposing that each one of the tetruples 
[ei], [e 3 ], ipi([e 3 ]), ipii([e 3 ]) and [ej, [e 2 ], ^i(N), ^n([e 2 ]) identifies a single 
circumference. No further condition comes from ip 3 and by symmetry we 
obtain an equivalent condition with respect to ip2- 

When s = another direct calculation shows that if u < 1/5 then t\ 
is not anymore the eigenvector of f\ corresponding to its largest eigenvalue. 
When u = a the circles r 13 and T 12 are tangent to each other and for u > a 
they intersect inside T A . □ 

Example 12. Let us give a short survey of the kind of geometry we meet in 
case of complex projective Sierpinski gaskets symmetric with respect to the 
imaginary axes. For u = 16/25 ~ a we get the gasket 

. 1 / 20 12A 1 (A A\ 1/4-4 

Al -7m{-^ 29j' A2 -7Hl v l 7)> A3 -^m{-l 7 



In Fig. 4 we show the orbit of a point under the action of the semigroup 
A generated by the A*. The triangle Ta is convex and, correspondingly, 
the triangle Za — T a \ (U^ =1 TaaJ is concave. Each angle is almost zero 
because the sides of the triangle are almost tangent to each other, which 
corresponds to the fact that the limit value a is close to 16/25. The restriction 
to Ta of corresponding maps ipi are contractive so Theorem 5 applies. A 
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a 



b 




Figure 4: Limit sets of complex self-projective Sicrpinski gaskets: (a) u — 16/25, s A — 
2.88, dim R A ~ 1.44; (b) u = 1/2, s A = 21og 2 3, dimi? A = log 2 3; (c) u = 9/25, s A ~ 2.88, 
dimi?^ ~ 1.44; (d) u = 1/5, ~ 2.61, dimi?^ ~ 1.305. In figure for each case we show 
the 19683 points of the orbit of a random point under the action of all matrices Aj of the 
gasket with |/| = 9. 
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rough numerical evaluation of the exponent of A gives sa — 2.88, so that 
dimR A ~ 1.44. 

By increasing u the curvature of the sides increases ( we consider negative 
the curvature of concave sides) until it gets zero for u = 1/2. The semigroup 
is now generated by 

Al = 7f(o O'^TlG 2)'^ = 7! (0 2 1 )- 

In this case all sides are segments of straight lines and the gasket is diffeo- 
morphic to the standard Sierpinski gasket in IR 2 . It is easy to prove that sa = 
21og 2 3 and, correspondingly, we get the well-known result dimi?^ = log 2 3. 

By increasing u further the curvature of the sides keeps increasing and 
therefore T4 becomes convex. For u = 9/25 the semigroup is generated by 

1 /3 6i\ . 1 / 3 3\ , 1 (Z -3\ 

T/ie corresponding ipi are contractive over T4 so £/ia£ Theorem 5 still ap- 
plies. A rough numerical evaluation of the exponent gives sa — 2.88 so that 
dimi? c ~ 1.44. 

At the extremal value u — 1/5 every angle of the triangle is equal to n, 
namely every triangle Za is actually a circle. Indeed this Sierpinski gasket is 
actually the Apollonian gasket A 3 , introduced in the Motivational Example 2 
and generated by 

T/izs tzme i/ie corresponding ipi are only non- expansive, which corresponds to 
the fact that A 3 are parabolic. 

Finally we point out that all these gaskets are fast. Here we outline the 
argument in case of the Apollonian gasket A 3 but the same argument holds 
for all complex projective gaskets symmetric with respect to the imaginary 
axes. Note first of all that it is straightforward proving by induction that 
|| Ai || = I 22 1 f or ever V matrix Aj e A 3 . Now consider the case 
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Then 

1 / a + c b + d 
23L ~ 2 \a + 5c b + 5d 

so that \\A I23L \\ = l\-f(b + d) + 5(b + 5d)\ > 2\6\\d\ > l\\A I \\\\A 23L \\. The case 
of A 32 is completely analogous to this. The remaining four combinations are 
instead analogous to the case of 

a 1 ( ~ * ^i 
12 = 2\-2 + i Q + i 

This time 

. If —ia + 3ic —ib + 3id 

12L = 2 V(-2 + i)a + (6 + i)c (-2 + i)b + (6 + i)d 

and \\A I12L \\ = ±\ 1 (3id-ib) + 5((i-2)b + (6 + i)d)\ > 2\5\\d\ > ||A 12L | 
Hence A 3 is a fast gasket with coefficient not smaller than 1/5. 



4.3 n>3,K = R 

In the real case, projective maps induced by at least 3x3 matrices are not 
conformal and we could find any simple way to prove analogues of Theorems 
4 and 5. The non-triviality of the matter is granted by the well-known non- 
triviality of the theory of real self-affine sets (e.g. see [Fal88, FL98, ABVW10, 
FM11]). Indeed, since PS'L^(M) contains a subgroup homeomorphic to the 
(n— l)-dimensional affme group, self-projective sets are at least as non-trivial 
as the self-affine ones (see Section 4.3.1 for more details). 

Because of this and in order to provide motivation for the interest of real 
self-projective sets we restrict our attention to the following particular case: 

Definition 6. Let F = (f±, . . . , f n ) be a free semigroup of volume-preserving 
linear automorphisms ofW 1 and ipi, . . . , ip n G PSL n (M.) the induced projec- 
tive automorphisms ofM.P^ 1 . Let S = {e±, . . . ,e n } be a n-frame of W 1 
and S* = {ei, . . . ,e n } its dual frame. We say that F is a real projective 
Sierpinski gasket over £ if the following conditions are satisfied: 

1. fi = A%e k <8> e? with A% > 0; 

2. /(e^ = Xfii, with \ = maxi< j <„{yl^}; 



44 



3. fi(ej) = ati + (3ej with a, (3 > 0; 

4- ^MM) = ^j(N), i^j- 

We say that S is a proper frame for F . More generally given m < n of 
the fi we say that they are a Sierpinski gasket if there exist automorphisms 
frn+i, ■ ■ ■ , fn such that (f±, . . . , f n ) is a Sierpinski gasket. 

Note that conditions 1-3 above imply that span{ei} is the only eigenspace 
of fi corresponding to its largest eigenvector, so that every proper frame for 
F identifies the same n points [e*] on RP n_1 . 

Denote by C(£) the positive cone over c, namely the convex hull of the 
set U™ =1 {Aej,A > 0}. Then its projection on MP™" 1 is the same for every 
proper frame of F and we denote it by Tp. This set is a (n — l)-simplex 
with the n points [e«] as vertices. By points 2 and 3 of the definition above, 

[e«] is a fixed point for ipi and each set Tpf^ip^Tp) is a (n — l)-simplex 
having in common with every other Tj-/., i ^ j, the vertex ^(fc?'])- Like 
in case of the (n — 1) -dimensional standard Sierpinski gasket in IR™ -1 , the 
difference between Tp and \Jf =1 Tpf. is the interior of a convex polyhedron 
with n(n — l)/2 vertices that we denote by Zp. 

By repeating this procedure recursively we see that, at every step k > 0, 



It is standard to call Rp = P\k>oTk,F the residual set of F. 

For sake of simplicity and conciseness we limit our discussion to the fol- 
lowing subclass of Sierpinski gaskets: 

Definition 7. We say that a Sierpinski gasket F = (f\,...,f m ) is sim- 
ple when each fi either has only one eigenvalue (first kind) or has exactly 
two eigenvalues and the eigenspace corresponding to the larger one is 1- 
dimensional (second kind). 

Example 13. The most important 1 -parameter family of real projective Sier- 
pinski gaskets we discuss in this paper is F" = {/",...,/"}, a > 1, 



T k /^{jT Ffl =Tp\ (JS 



|/|=fc |7|<fc 
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For n = 3 the f" are represented, with respect to any proper frame, by 
the matrices 

fa 1 1\ (\ 0\ /l N 

A?' 3 = a-3 1 Uf = «-Hl « 1 Lf = «4 1 
\0 1/ \0 1/ \1 1 a, 

yls already shown in the introduction, for a = 1 we get the cubic gasket C 3 
and for a = 2 the (real projective generalization of the) standard Sierpinski 
gasket S 3 . 

Consider now the dual semigroup F* C Aut((M, n )*) of a Sierpinski gasket. 

Proposition 11. Let F be a simple Sierpinski gasket over a n- frame £ = 
{e±, . . . , e n } of M. n generated by maps 

Then F* is a simple Sierpinski gasket over the frame H = {i] 1 , . . . ,i] n } of 
(R n )* defined by 

v i = j2^ £j + ^-^y i - ( 25 ) 

Proof. By direct calculation we see that 

f:W) = £%/^ J ) + («i - 7i)7?(e*) = 

= ^ + ~ 7*)(ai£ < + ^ hi 6 ?) = 



and 



f:(v k ) = £/v^') + ( a k-ik)f*(£ k ) 



3+k 



□ 
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Lemma 3. Let F = (/,) C Aut{V n ) be a simple Sierpinski gasket over 
£ = {ei}. Then the following inequalities hold: 

\\fi\\fi < C min {||//(e fc )|^i + ||//(e fc /)||^i}, if F is of the first kind. 

l<k,k'<n 
kj=k' 

\\fi\U 1 < C min \\fi(ek) ||^i, if F is of the second kind, 

l<k<n 

(26) 

for some C > 0, where \\uj\le 1 = J2i<j< n \ u j\- 

Proof. Let H = {rf = ^2$jE^} be the proper frame for F* introduced in 
(25), where /3j = > 0, j ^ i, and f}\ = — 7i > 0. Clearly u = 
^2 1<i<n rf e C{%). By the previous proposition, uij = fj(u)) G C("H) for all 
I e l n . This means that ui = Y.i<i< n ( u i)i e% = Z)i<i<n'W ( witn \ > 
and J2i<i< n > °) so that = Si<j<„ and therefore 

^~ ^ _^ ^ 

l<j<n l<*ii<i l<i<n 

Now note that (a>j)fc > (min^g (3\) S/?|>o ^* ls & l wa y s a non-empty con- 
dition. 

If F is of the second kind then /3j > for all i, j, so that 

max {/3H 

ll n , l<»J<n / x 

< n — , a K) fc 

l<j<n 

for all 1 < k < n. 

If F is of the first kind then ft l k = iff % — k, namely every w fc misses 
Afc in its expression, which we can recover by adding any other Uk>, k' ^ k. 
Hence in this case 

max l A 

\\m\\ <{n- i) " '," 3 . (Mk + (ui)k') 

mm {/?]} 

for all 1 < fc, fc' < n, k ^ k! . 
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Finally note that // = Xa<jj<n^/j e * ® £ * ■> so that 

l<i,j,fc<n 

and therefore 

M k < max E Aj fc = max {^.}||/?(e fc )||^ 

and 

(min^OHMI^ = (min/^E^ < IM* 

from which follows the claim of this lemma. □ 

Example 14. Consider the Sierpinski gaskets F a ' n introduced in Exam- 
ple 13. A proper frame for (F Q,n )* is given by rf = (a — \)e l + Ylj^i^ > 
so that 00 = Y^i<i< n r t — (a + n — 2) Yl<i<i<n e% and therefore 

W , = /?( W ) = (a + n-2) J] /?(e*) = (a + n-2) £ ^ fc = 

l<i<n l<j,fc<n 

= (a + n-2) J2 Wfi(e k )\\^ k . 

l<k<n 

If a = 1 t/jen 

i,3 3 

and 

Mi = A ^ 

so that 

< (n - 1) ((w/) fc + (a;/) fc /) 

or, equivalently , 

||//||/i<(n-l)(||//(e fc )|| /1 + WfMWe) 

for every k ^ k' 
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If a > 1 then 



\H = (n - 1) max{a - 1, 1} A, 




and 



> minjo; 



1.1} E A i 



SO t/iflt 



^/ll^ 1 < — 1) min{a; — 1, 



or, equivalently, 



\\fj\y < (n-l)irun{a-l,l}||/j(e fc )||*i, 



for all k. 

Proposition 12. Every simple Sierpinski gasket is a fast gasket. 

Proof. The key fact here is that in every Sierpinski gasket F with a proper 
frame S, for any i ^ j and every k, fij(e k ) is linearly dependent on both 
and ej. Indeed if k ^ j then 

fij(e k ) = fi(Pjkej + ajkGk) = ot jk fi(e k ) + f3 jk ((5^ + a^ej) , 

while if k = j then 



Hence the matrix representing with respect to £ has at least (actually, 
exactly) two rows with all non-zero coefficients. This means that every col- 
umn of the matrix representing fijL is a linear combination with strictly 
positive coefficients of /z,(ej) and fh{^j) with possibly some other positive 
contribution from the other vectors. 

From this we deduce immediately that 



for some C > 0. Since > C'll/yxll^ 1 for all i ^ j and some C" > 0, 



W/iulU >C\\fi\W (11/^)11* + 11/^)11,!) 



for some C > and therefore, by Lemma 3, that 



\\fii jL \W >c"||MHIMk 



our claim follows. 



□ 
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Proposition 13. Let F C AutiW 1 ) be a simple Sierpinski gasket, £ a proper 
frame of F and /x any measure ofM.P"^ 1 in the measure class of the round 
measure, namely the measure induced by the metric of sectional curvature 
identically equal to 1. Then there exist constants A, B,C, D > s.t. 

< KTr fl ) < Trrinr-, TO < ^ Ffl ) < 



\\fi\\ n -^ 111 - iimi-' \\fi\\ n -^ 111 - \\h\r 

where a n = n if F is of the second kind and a n = n — 1 if it is of the first 
kind. 

Proof. It is enough to prove the claim in some chart containing Tp. We fix 
coordinates (x 1 , . . . , x n ) so that the vectors of S are 

e 1 = (l,0,...,0,l),...,e n _ 1 = (0,...,0,l,l),e n = (0,...,0,l) 

and use the chart x n = 1. In this chart we pick any smooth measure v of 
finite total volume and with constant density equal to 1 within Tp- 

Note that // = A\ k Ci®e k = A l Ik e{dj®e k , where the last row of the matrix 
A\ k e\ contain the i 1 norms of the vectors //(ej). A direct calculation shows 
that 

1 1 

K T Ff,) - n - n • 



ni 

k=l k=l 



II WfiMWe 

Clearly k ~ 1 < 1. By Proposition 3, if F is of the second kind then 



mi/. 



i{e k ) \\ e i 



A < — 



IIMI-i 

for some A > 0. If it is of the first kind assume, for the argument sake, that 
||//(e 1 )||,i<...||//(e n )||. Then 



< 



nn/'( e * n n//(e*)ii^ 

k=l k=l 
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<c TT ILMgfc)]^ + \\fi(e k +i)\\ii < 2 „_ lc 

The geometry of is more complex. We divide it in n {n — l)-simplices 
Zj, where Z^s vertices are the (n — 1) points [fi(fi(ej))], j 7^ h plus the point 
Ki<i<„ //(*)]■ Then 

1 1 



l<i<n jy^i 

since all components of all vectors are positive. Even in this case then 

1 

< n\ 



»{Zi)\\fi\\l 
and 



\\fi\\pY\rtiWij{ e i)i + li( e i)j\\p ~ max {/ 3 >7}n i ^ll//(e J ) + /(e J )l|£i 

for some A. □ 

Next proposition supports the idea that the residual set of a Sierpinski 
gasket, or at least of a simple one, have non-integer dimension: 



Proposition 14. The residual set of a simple Sierpinski gasket F C SL^ , 

yn— 1 



-± 

is a null set with respect to the measure class of the round measure on MP 
Proof. From the previous proposition we see that 



1 < KT Ffl ) < c 

K Z F fl ) \\fl I 



if F is of the first kind and 



1<^4<C 



if it is of the second. Let us first assume that F is of the second kind and let 

S k = J2 M^/J and P k = J2 ^Z Ffl ). 
\i\=k+i \i\=k 
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Then 

Sk+i = Sk — Pk < Sk(l — C) 

and therefore 

Sk — Sk+i > CSk 
so that, after making a telescopic sum, we get 

Si — lim Sk > C lim /cS^ 

which immediately implies that lim^oo Sk = 0. 

If F is of the first kind then /x(7jp/ / )///(Z F / / ) < C/ min|j| = |/| ||/j||. It is 
easy to check that min|j|=|/| is proportional to |/| and therefore 

S k+ i = S k -P k < S k (l - C/k) 

and 

Sk — Sk+i > CSk/k 

so that 

5i - lim S k > C lim S k V 1/j. 
Since the series 1/j diverges we get again that lirn^oo Sk = 0. □ 



4.3.1 AfRne Sierpinski Gaskets 

In [FL98] Falconer and Lammering studied in detail the family of affine 
Sierpinski gaskets S a> b, a,b e (0, 1), defined by the affine transformations 




Notice that the case a = \,b — | corresponds to the standard Sierpinski 
gasket. In particular they proved that the box dimension of the corresponding 
residual set R a ^ is given by the unique root of the equation 

(1 -ay + ab"- 1 +a(l - h) 3 ' 1 = 1 (27) 
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in the triangle 7\ = {(a, b) G (0, l) 2 , a > max{6, 1 — b}} and of 

(1 - a) s + a 8 ' 1 = 1 (28) 

in the opposite triangle T 2 = {(a, 6) G (0, l) 2 , a < min{6, 1 — 6}}. This setting 
provides a convenient source of examples for comparing the box dimension 
of the residual set of a real projective Sierpinski gasket with the correspond- 
ing gasket exponent. Indeed the injection sending the affine transformation 
S(x) = Tx + v, with T G M 2 (R) and x,u G M 2 , into the 3 x 3 matrices 

S — > ^ ^ applied to the gives the three matrices 

'1-6 l-a-6 o N 

Uj = | I - a a 1 . M 2 = | a | . VA; = | a 

1, 

so that the action induced by the on MP 2 in the affine chart z = 1 coincides 
with the action of the Si. 

Like in Example 9, these matrices are upper triangular and ||M/|| = 1 for 

(a X /A 
/3 v \ and notice that 
lj 

< A + (3 < 1 and < // + v < 1. Indeed we can limit the discussion to 
left multiplication by Mi and M 3 and it is immediate to verify by induction 
that, assuming the inequalities above for Mj, the products 

((1 — a)a (1 — a)X (1 — a)fj, 

(l-a)/3 (l-a)u + a 

1 

and 

/(l-b)a (l-b)X + (l-a-b)/3 (l-o)/i + (l-a-o)z/ + 6 N 

M 3 M = a/3 az/ 

\ 1 

satisfy the same inequalities. Now consider the semigroup A a b generated by 
the At = Mi/detM- /3 G SL 3 (R). Clearly = det M" 1 ^ 3 • • • det M~^ 3 

for every I — %\ ■ • • and therefore 

U a , b ,k = £ ll^ll" 5 = ( det M l /3 + det M 2 /3 + det M 3 /3 )"- 
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Figure 5: Affine Sierpinski gaskets A a ^ for four possible pairs a, b. For each 
one we plot (in green) the set T 7 ^ ab . For the upper two the box dimen- 
sions can be evaluated analytically and their first six digits are, respectively, 
dime -R 3 1 = 1.72368 and dim^i?i i = 1.68886. The numerical evalua- 
tion of the box dimensions with an elementary box-counting algorithm gives: 
dime i? 3 i ~ 1.71, dime i?i i — 1.66, dim^-Ri i — 1.60, dim^ i?i 3 ~ 1.60. 

4'2 4 ' 2 2 >4 2 '4 

See Table 1 for a comparison of the box dimension of these and other affine 
gaskets with the expoentn of the corresponding real projective Sierpinski 
gaskets. 
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By Theorem 2 then the exponent SA a b is the unique solution of the equation 

(1 - a) 2s ' z + (ab) s / 3 + (a(l - 6)) s/3 = 1. (29) 

Proposition 15. 2sA ab < 3(Wm B R ab within the two triangles T li2 , £/ie 
egna/ sign holding only in their common vertex. 

Proof. After writing (29) in terms of t = 2s/ 3 and renaming t to s we are 
left with the equation (1 — a) s + (ab) s ^ 2 + (a(l — b)) s / 2 = 1. Comparing this 
expression with the left-hand sides of (27) and (28) we see that it is enough 
to prove that 

(ab) s ' 2 + (a(l - b)) 3 ' 2 < minK- 1 , ab^ 1 + a(l - by- 1 }. 

Since for obvious geometrical reasons dims R a ,b < 2 we can assume in the 
following s > 2. Let us denote respectively by f a ,b{s) , g a (s) , h a ^{s) the three 
functions above and notice that, since by hypothesis a, b, 1 — b G (0, 1), they 
are all strictly monotonically decreasing functions of s converging to as 
s — )• oo. Moreover f a ,b(2) = <7a(2) = h a> f,(2) = a and since these functions can 
have only one intersection it is enough to verify their behaviour for s — > 0. 
A direct calculation shows that, for every pair a,b e (0, l) 2 , we have that 
lim s _> fa,b(s) = 2 while lim s ^ g a (s) = lim s ^ h a (s) = oo. □ 

Numerical experiments (see Fig. 5) clearly suggest that, even outside of 
the triangles Ti 2 , dims R a ^ always larger than \sA a b with the only exception 
of the case a = 1/2, b = 1/2, when these two quantities coincide. Moreover 
it appears that, roughly, \sA ab > jjj dime Ra- 

4.3.2 The cubic semigroups C" 

Recall that, by definition, /"(e*) = cce^, /"(e^) = + e^, j ^ i. 

n = 3. In IR 3 we use coordinates (x, y, z) with respect to the frame e[ = 
ei + e3, e' 2 = e 2 + e^, e' 3 = e 3 , so that the fi are represented by the matrices 

(a - 1 1\ /l 0\ / 1 0\ 

Ai= 1 , A 2 = a - 1 1 , A 3 = 1 . 
\a- 2 2/ \0 a-2 2/ \2-a 2-a a) 

In the affine chart [x : y : z\ — >■ (w, t> ) = (x/z,y/z) of IRP 2 we have therefore 
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a 


b 


S A 

° A 

(num.) 


S A 

(anal.) 


2s A /3 


dim d /? a 
(num.) 


dime /? a 
(anal.) 


1/4 


1/2 


2.44 


2.42632 


1.61755 


1.66 


1.68886 


3/4 


1/2 


2.48 


2.45425 


1.63617 


1.71 


1.72368 


1/2 


3/4 


2.35 


2.34443 


1.56295 


1.60 




1/5 


3/10 


2.44 


2.43735 


1.62490 


1.76 


1.71262 


4/5 


3/10 


2.47 


2.46960 


1.64640 


1.77 




3/10 


1/5 


2.43 


2.37354 


1.58236 


1.75 




7/10 


1/5 


2.35 


2.35249 


1.56833 


1.72 


1.63373 



Table 1: Values of the exponents of affine Sierpinski gaskets A a ^ for several 
pairs a, b and of the box dimension of the corresponding residual sets R a ,b- 
Numerical evaluations for sa were done by calculating N A (k) for the values 
k = 2 P , p = 1, . . . , 12, and are presented to motivate our confidence in a rel- 
ative error not bigger than 1% in the other evaluations provided throughout 
the paper when an analytical evaluation is not available. Numerical evalua- 
tions for the box dimension of R a ^ were done via an elementary box-counting 
algorithm and a comparison with the available analytical evaluations suggest 
that their relative error is about 10%. 



that 



(j; ( ) = ( (a-i)u+i v \ 

flK^iil) \ja-2)u+2' (a-2)u+2j 



-a) 



x)(u+v)+a'> (2— a)(u+v)+c 

and the vertices of the invariant triangle T F a are [ei] = (1,0), [e 2 ] = (0,1) 
and [e 3 ] = (0,0). A direct calculation of the eigenvalues of the Jacobian 
matrices Dipi shows that, within Tp«, 

1 (X 1 Oi 

min {-> -r)d{x,y) < d(ipi(x)^i(y)) < max{-, -}d(x,y) 
a A a A 

for all i = 1,2,3, namely the semigroup (ipi,^,^), as a IFS, is hyperbolic 
for a G (1,4) and parabolic for a = 1,4. The ipi are not contractions with 
respect to the Euclidean distance in this chart for the other values of a (see 
Fig. 6 for the plot of T^^a for several values of a). 

Analytical bounds for the Hausdorff dimension of the residual sets 
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a 


Si?" 


2sf«/3 


dime Rf" 


1 


2.447 


1.631 


1.72 


1.3 


2.395 


1.596 


1.72 


1.7 


2.377 


1.585 


1.71 


2 


2.359 


1.573 


1.59 


3 


2.378 


1.586 


1.71 


4 


2.389 


1.593 


1.73 


7 


2.394 


1.596 


1.76 



Table 2: Numerical evaluation of the exponent of the real projective gaskets 
F£ and of the box dimension of the corresponding residual sets for several 
values of a. No analytical formula is known for these quantity. These data 
confirms the relation 2s^/3 < dim^-R^i already observed in Table 1 and the 
fact that roughly 25^/3 > 9dim B i?^/10. 

can be obtained via Propositions 9.6 and 9.7 in [Fal90], namely 

4 1 4 1 

min{log 3 -, } < dim H R F? < max{log 3 -, }. 

a log 3 a 3 a log 3 a 

For a = 2 we get, as expected, dim# Rp» = log 3 2. 

Analytical bounds for the exponents sf« can be obtained from Theorem 1. 
Here we present calculations for F 3 = C 3 , the cubic gasket. Due to the 
symmetry between the generators it turns out that 

Hc 3 (s) = 6/i C3 A 12 = 3 • 2 1 - s ((s), 

from which, as the unique solution of Hc :i { s ) = 3 s , we get the lower bound 
1.52 < sc 3 . To get the first upper bound we must consider the function 

KM = 3 • 2^ fa) + 2 2 - s ((s, ? -) - 2 1 - - 3^ , 

from which we get 1.7 < sc 3 < 7.1 as the unique solutions of //c 3 ,2(s) — 3 ±s . 
In order to get more meaningful bounds we should consider some /J,c 3 ,k with a 
large k but leave this to a future paper. Interpolating on the curve log Nc 3 (k) 
as function of log A; for k = 2 r , 1 < r < 13, we get a reliable estimate of 
sc 3 c± 2.444. A rough numerical evaluation of the box dimension of R Ci by 
counting the number of squares needed to cover the fractal gives dim# ~ 1.72, 
compatible with the relation 3 dims Rc 3 > 2sc 3 suggested in Conjecture 1. 
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7 10 100 

Figure 6: Real projective Sierpinski gaskets for several values of a. For 
each one we plot (in green) the set T 7 ^a. Heuristic numerical estimates of 
their exponents and of the box dimension for the corresponding residual sets 
for a < 7 are given in Table 2 
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n > 4. In R n we use coordinates (x 1 , . . . , x n ) with respect to the frame e' x = 
e„_i + e n , e' n = e n . For n = 4 the matrices A\ and A 4 are 



e\ + e ri 
given by 



fa -I 1\ 

10 

10 

\a-2 2/ 



A 4 



/ 1 0\ 

10 

10 

\2 — a 2 — a 2 — a a J 



and A 2 and A 3 can be obtained via permutations of A 1 . Similarly happens 
for n > 4. Correspondingly we use coordinates u l = x l /x n , % = 1, . . . , n — 1, 
and obtain 



1W + 2 



it 



n-l 



,a 



l)^ 1 + 2 ' (a - IV + 2 ' " ' ' (a - l)^ 1 + 2 
similarly for i < n — 1 and 



^n-l(^) 



n-l 



(2-a)(w 1 H hw"" 1 ) + (2 -a){u l H htt"" 1 ) +a 



A direct evaluation of the eigenvalues of the Jacobian matrices of the ipi gives 
the same result we got for n = 3. In particular for every n > 3 we have that 
the gasket F% is a hyperbolic IFS for a G (1,4) and a parabolic IFS for 
a = 1,4. The bounds on the Hausdorff dimension of the residual sets give 



min {!og„ 



a ' log n a 



1 4 
} < dim// Rpa < max{log n 



a ' log n a 



}• 



For n = 4, ct = 2 we get the well-known result that the dimension of the 
standard Sierpinski tetrahedron is equal to 2. Numerical evaluations suggest 
that the same could hold for the 4-dimensional version of the cubic gasket 
(see Fig. 7 for a picture of the two sets). 



4.3.3 The Apollonian gasket 

We conclude the paper with a brief discussion on the Apollonian semigroup, 
namely the semigroup H C SX 4 (N) generated by the matrices H 1 ,H 2 ,H 3 
introduced in the Motivational Example 2. This case was thoroughly stud- 
ied, somehow implicitly, by Boyd, in particular in [Boy 72, Boy73b, Boy82], 
in the context of the sequence of curvatures in an Apollonian gasket and 
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Boyd's investigation and arguments were the archetype for most results and 
arguments in Section 3 of the present paper. 

Recall that the matrix Hi has only one eigenvalue (hence equal to 1) and 
therefore, even though H 2 and if 3 have eigenvalues larger than 1, H is a 
parabolic gasket. Next proposition grants that our results do apply in fact, 
as expected, to H itself and shows that s H < oo with arguments that entail 
only the matrices themselves. 

Lemma 4. Assume that matrices A 1 ,...,A m e M n (N) have the following 
properties: 

1. they have some number k ^ 1 of rows containing a single entry equal 
to 1 and all other equal to and these entries equal to 1 belong all to 
different columns and in those columns all entries are either or 1; 

2. other rows are such that each of their entries is smaller than the sum 
of the remaining n — 1 entries. 

Then this property is shared by all products of the Ai. 

Proof. We prove the lemma by induction. It is enough to consider the prod- 
ucts of two generic matrices A = (Aj), B = (Bj), satisfying the hypotheses. 
Assume first that k = for B, namely J2k^e B l — B \ for all i,£. Then 

B^)i= E m = E 4I>£> E m = ( AB Ye- 

k±t l<j<n l<j<n kj^l l<j<n 

k+l 

Assume now that k > 1 for B and denote by / = (ii, . . . ,ik) the rows 
with a 1 and all other entries equal to 0. Every line (if any) of A with a 1 
and all other entries equal to leaves unaltered the corresponding row in B 
and therefore the new line satisfies the conditions in the theorem. Otherwise 
we notice that 

J2(AB)l= £ A)B{> £ A)B> + 4E^" 

k+t l<j<n l<j<n l<j<n k^l 

k+l j&i jei 

If j e I then J2k^e B l * s either or 1. Since by hypothesis there are at least 
two such rows and Y^k^e B l + Y.kM B ( ^ 1 for eveT Y 3 J' e ^ 3 ^ /> and 
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the corresponding entries A 1 - and A 1 -, are both equal to 1, then 
E A )J2 B l^ E A)Bi = (AB)\, 

l<j<n k^i l<j<n 

jei jei 
therefore T. k ^B)\ > {AB)\. □ 

Proposition 16. The Apollonian semigroup is a fast gasket with coefficient 
c > 1/4. 

Proof. Note first of all that Hirst matrices satisfy previous Lemma's condi- 
tions. Moreover the entries in the third line are not smaller than all other 
entries in the same column and it is easy to see by induction that this prop- 
erty is preserved by products. 

Let Halloo = maxi<j< n Yli<j< n l^jl- A 1°°^ & t the 6 matrices Hij, i ^ j, 
shows that their third column has always at least three non-zero entries, so 
that ||Atj||oo > ||A/||oo Hj^j \^%\ where jo is the index of the element of 
the third column (if any) equal to zero (otherwise just set j — 1). By the 
previous Lemma and the fact that the norm of every A e H in concentrated 
in the third row, the sum of any three entries of the third row of A is always 
larger than || A ||, so that j || oo > ||A/||oo||^j||- Since 4||^|| > WA]]^ > \\A\\, 
the claim follows. □ 

Analytical bounds for the exponent sh of the Apollonian semigroup were 
studied in detail by Boyd in [Boy 70, Boy 72, Boy 73a] and we do not attempt 
to improve them here. 

Increasingly accurate numerical evaluations of sh with several different 
techniques have been given over the last half-century by Melzak [Mel69], 
Boyd [Boy82], Manna and Herrmann [MH91], Thomas and Dhar [TD94] and 
McMullen [McM98] giving respectively the following values, with a heuris- 
tic error of 1 unit on the last digit: 1.306951, 1.3056, 1.30568, 1.30568673, 
1.305688. We remark that, among all these evaluations, the one with the 
largest number of digits, given by Thomas and Dhar, is the only one based 
on a heuristic method, while the others are based on exact methods. 

Partly to test our own software evaluating the function Nn{k) for a 
generic gasket H and partly because the computational power of comput- 
ers increased quite a lot over the last fifteen years, which is how old is the 
last evaluation of the exponent, we repeated the elementary evaluation made 
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c 2 

(19) 


3, 15, 71, 287, 1231, 4911, 19831, 79279, 318383, 1273807, 5098247, 
20391887, 81590055, 326364583, 1305483999, 5221928631, 
20888160751, 83552534287, 334211194663 


c 3 

(13) 


4, 22, 148, 760, 4594, 24646, 136372, 740650, 4046188, 22022770, 
119929126, 652445212, 3550689778 


c 4 

(12) 


5, 37, 293, 2197, 15125, 103669, 714245, 4849045, 32901077, 
222724789,1507986917, 10202765749 


A 3 

(13) 


3, 12, 64, 316, 1784, 10004, 58224, 341386, 2033906, 12170708, 
73208110, 441772966, 267292497 


H 

(39) 


0, 1, 3, 8, 18, 48, 113, 278, 681, 1722, 4238, 10488, 25927, 64086, 158266, 391062, 967315, 
2390800, 5909752, 14608522, 36115118, 89275994, 220684802, 545546400, 1348603780, 
3333755028, 8241076212, 20372155276, 50360227721, 124491161884, 307744098990, 
760747405278, 1880578271904, 4648814463680, 11491932849933, 28408221038996, 
70225503797745, 173598409768852, 429137646728801 


(35) 


2, 7, 16, 34, 84, 151, 348, 679, 1546, 3034, 6546, 13476, 28409, 59578, 122139, 261698, 
531191, 1144823, 2314772, 4986951, 10132768, 21667197, 44400099, 94074745, 194587388, 
408651488, 852101402, 1777247239, 3726410796, 7738675037, 16274400897, 33739772516, 
71002774691, 147235829060, 309533001058 



Table 3: Values of Na^) for small k for the cubic semigroups Ci, % = 2, 3, 4, 
the Apollonian semigroup A3, the Hirst semigroup H and the semigroup F 
of Example 10. In the left column it is also reported the numnber of terms 
displayed in the right one. 



by Boyd in 1982 by evaluating N H {k) for k = 2 P , p = 1, . . . , 40 with re- 
spect to the norm \\Aj\\ = Xa<« j<i{A-i)ijV l w l , where v = (—1,2,2,0) and 
w = (1,1,1,2) (this way \\Aj\\ is equal to the the curvature of the circle 
of multi-index / in the Apollonian gasket generated by the circles of radius 
— 1,2,2), and then interpolating the data obtained (see Table 3). We found 
a value of sh — 1.30568673 which fully confirms the heuristic evaluation of 
Thomas and Dhar and suggests an error of 2 on the last digit of the estimate 
of McMullen. 

Acknowledgments 

I am grateful to S.P. Novikov, LA. Dynnikov and B. Hunt for several pre- 
cious insights and suggestions on the matter presented in this article and 
thank them and T. Gramchev for several fruitful discussions while writing 
the paper. Finally I am grateful to my wife M. Camba for helping speeding 
up considerably my code to evaluate the functions N(k). Most calculations 
were done on the ~100 cores 2,66GHz Intel Xeon Linux cluster of INFN 



63 



(Cagliari); latest calculations were also performed on the iMac cluster of 
the Laboratory of geometrical methods in mathematical physics (Moscow), 
recently created by the Russian Government (grant no. 2010-220-01-077) 
for attracting leading scientists to Russian professional education institutes. 
Finally I an grateful to the IPST and the Mathematics Department of the 
University of Maryland for their hospitality in the Spring and Fall 2011 while 
I was working at the paper. 

References 

[ABVW10] R. Atkins, M.F. Barnsley, A. Vince, and D.C. Wilson, A charac- 
terization of hyperbolic affine iterated function systems. Topology 
Proc. 36 (2010), 189-211. 

[Boy70] D.W. Boyd, Lower bounds for the disk packing constant. Math. 
Comp. 24 (1970), 697-704. 

[Boy71] , On the exponent of an osculatory packing, Can. J. Math. 

2 (1971), 355-363. 

[Boy72] , The disk packing constant, Aeq. Math. 7 (1972), 182- 

193. 

[Boy73a] , Improved bounds for the disk packing constant, Aeq. 

Math. 9 (1973), 99-106. 

[Boy73b] , The residual set dimension of the Apollonian packing, 

Mathematika 20 (1973), 170-174. 

[Boy82] , The sequence of radii od the Apollonian packing, Math- 
ematics of Computation 39 (1982), no. 159, 249-254. 

[BV10] M. Barnsley and A. Vince, Real projective iterated function sys- 
tems, Journal of Geometric Analysis (2010), 1-36. 

[Cox37] H.S.M. Coxeter, The problem of Apollonius, The American 
Mathematical Monthly 75 (1937), no. 1, 5-15. 

[DD09] R. DeLeo and LA. Dynnikov, Geometry of plane sections of the 
infinite regular skew polyhedron {4, 6|4}, Geometriae Dedicata 
138 (2009), no. 1, 51-67. 



64 



[Fal88] K.J. Falconer, The Hausdorff dimension of self-affine fractals, 
Math. Proc. Camb. Phil. Soc. 103 (1988), 339-350. 

[Fal90] , Fractal geometry, Wiley & Sons, 1990. 

[FK60] H. Furstenberg and H. Kesten, Products of random matrices, The 
Annals of Mathematical Statistics 31 (1960), no. 2, 457-469. 

[FL98] K.J. Falconer and B. Lammering, Fractal properties of general- 
ized Sierpinski triangles, Fractals 6 (1998), no. 1, 31-41. 

[FM11] K.J. Falconer and J.J. Miao, Local dimensions of measures on 
self-affine sets, 2011, arXiv:l 105.241 lvl [math.MG]. 

[GLM+03] R.L. Graham, J.C. Lagarias, C.L. Mallows, A.R. Wilks, and C.H. 

Yan, Apollonian circle packings: number theory, Journal of Num- 
ber Theory 100 (2003), no. 1, 145. 

[GLM + 05] , Apollonian circle packings: Geometry and group theory 

I. the Apollonian group, Discrete & Computational Geometry 34 
(2005), no. 4, 547-585. 

[GLM + 06] , Apollonian circle packings: Geometry and group the- 
ory ii. super- apollonian group and integral packings, Discrete & 
Computational Geometry 35 (2006), no. 1, 1-36. 

[Hir67] K.E. Hirst, The Apollonian packing of circles, J. of London Math. 
Soc. 42 (1967), 281-291. 

[Kir07] A. A. Kirillov, A tale on two fractals, 2007, math- 
proxy.sas.upenn.edu/~kirillov/MATH480-F07/tf.pdf. 

[KOll] A. Kontorovich and H. Oh, Apollonian circle packings and closed 
horospheres on hyperbolic 3 manifolds, J. of Am. Math. Soc. 24 
(2011), 603-648. 

[Lev93] G. Levitt, La dynamique des pseudogroupes de rotations, Inven- 
tiones Mathematica3ll3 (1993), 633-670. 

[McM98] C.T. McMullen, Hausdorff dimension and conformal dynamics, 
III: Computation of dimension, American Journal of Mathemat- 
ics 120 (1998), no. 4, 691-721. 



65 



[Mel69] Z.A. Melzak, On the solid-packing constant for circles, Math. 
Comp. 23 (1969), 169-172. 

[MH91] S.S. Manna and H.J. Hermann, Precise determination of the frac- 
tal dimensions of Apollonian packing and space-filling bearings, 
J. Phys. A: Math. Gen. 24 (1991), 481-490. 

[MU98] R.D. Mauldin and M. Urbanski, Dimension and measures for a 
curvilinear Sierpinski gasket or Apollonian packing, Advances in 
Mathematics 136 (1998), 26-38. 

[NovOO] S.P. Novikov, I. Classical and Modern Topology. II. Topological 
phenomena in real world physics., Modern Birkhauser Classics, 
vol. GAFA 2000, pp. 406-425, 2000, arXiv:math-ph/0004012. 

[01d96] A. Oldknow, The Euler-Gergonne-Soddy triangle of a triangle, 
The American Mathematical Monthly 103 (1996), no. 4, 319- 
329. 

[PA02] P. Pajares-Ayuela, Cosmatesque ornament: Flat polychrome geo- 
metric patterns in architecture, W. W. Norton & Company, 2002. 

[Pap40] Pappus of Alexandria, Collection, Book VII, ca. 340. 

[Sarll] P. Sarnak, Integral apollonian packings, American Mathematical 
Monthly 118 (2011), no. 4, 291-306. 

[Siel5] W. Sierpinski, Sur une courbe dont tout point est un point de 
ramification, C. R. Acad. Sci. Paris 160 (1915), 302-305. 

[Sod36] F. Soddy, The kiss precise, Nature 137 (1936), 1021. 

[Sul84] D. Sullivan, Entropy, haus dor ff measures old and new, and limit 
sets of geometrically finite kleinian groups, Acta Mathematica 
153 (1984), no. 1, 259-277. 

[TD94] P.B. Thomas and D. Dhar, The H aus dor ff dimension of the Apol- 
lonian packing of circles, Journal of Physics A: Mathematical and 
General 27 (1994), no. 7, 2257-2268. 

[Vinl2] A. Vince, Mobius iterated function systems, 2012, 
www.math.ufl.edu/~avince/Papers/MobiusAMS.pdf. 



66 



